Accepted for publication in the Astrophysical Journal 

Preprint typeset using I4TpX style emulateapj v. 1 1/10/09 



CONSTRAINTS ON THE ASSEMBLY AND DYNAMICS OF GALAXIES: I. DETAILED REST-FRAME OPTICAL 
MORPHOLOGIES ON KILOPARSEC-SCALE OF Z - 2 STAR-FORMING GALAXIES 1 

N. M. FORSTER SCHREIBER 2 , A. E. SHAPLEY 3 , D. K. ERB 45 , R. GENZEL 2 ' 6 , C. C. STEIDEL 7 , N. BOUCHE 4 ' 8 ' 9 , G. CRESCI 10 , R. 

Davies 2 

Accepted for publication in the Astrophysical Journal 



ABSTRACT 

We present deep and high-resolution Hubble Space Telescope NIC2 F160W imaging at 1.6 pm of six z ~ 2 
star-forming galaxies with existing near-infrared integral field spectroscopy from SINFONI at the Very Large 
Telescope. The unique combination of rest-frame optical imaging and nebular emission-line maps provides 
simultaneous insight into morphologies and dynamical properties. The overall rest-frame optical emission 
of the galaxies is characterized by shallow profiles in general (Sersic index n < 1), with median effective 
radii of R & ~ 5 kpc. The morphologies are significantly clumpy and irregular, which we quantify through a 
non-parametric morphological approach, estimating the Gini (G), Multiplicity (<£), and M20 coefficients. The 
estimated strength of the rest-frame optical emission lines in the F160W bandpass indicates that the observed 
structure is not dominated by the morphology of line-emitting gas, and must reflect the underlying stellar mass 
distribution of the galaxies. The sizes and structural parameters in the rest-frame optical continuum and Ha 
emission reveal no significant differences, suggesting similar global distributions of the on-going star formation 
and more evolved stellar population. While no strong correlations are observed between stellar population 
parameters and morphology within the NIC2/SINFONI sample itself, a consideration of the sample in the 
context of a broader range of z ~ 2 galaxy types (X-selected quiescent, AGN, and star-forming; 24 /im-selected 
dusty, infrared-luminous) indicates that these galaxies probe the high specific star formation rate and low stellar 
mass surface density part of the massive z ~ 2 galaxy population, with correspondingly large effective radii, 
low Sersic indices, low G, and high \ff and M20. The combined NIC2 and SINFONI dataset yields insights of 
unprecedented detail into the nature of mass accretion at high redshift. 



Subject headings: galaxies: evolution 
galaxies: structure 
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1. INTRODUCTION 

A fundamental but still unsolved issue is the origin of 
the Hubble sequence and the formation of the spheroidal 
and disk components of present-day galaxies. Observa- 
tions indicate that the Hubble sequence emerges around z ~ 
1-2, along with the closely connected bimodality in inte- 
grated rest-frame optical colors characterized by a "red se- 
quence" of mostly elliptical/bulge-dominated early-types and 
a "blue clo ud" of generally disk-dominated and irregular late- 
types (e.g. [Bell et al. 12004 Ivan den Bergh et al. 1996 J200H 
iLilly et all Il998t IStanford et al.l l2004T |Ravindranat h et al.l 
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20051 iKriek etal.l l2008bl 120091 
Although colors are easier to mea- 
sure, morphologies potentially yield more information about 
the processes that shape galaxies over time. Morphologi- 
cal studies reveal a progressively larger proportion of irreg- 
ular and disk-like types at hi gher redshift compared to the 
local galaxy population (e.g. jLotz. Primack, & Madau ll2004|: 
i Lotz et al. I liooa lLaw et all l2007at IScarlata et al. I I2007albl: 
Elme green et al. 1 120071: iPannella et al. 1 120091) , hinting at the 
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diversity of processes driving early galaxy evolution. 

Morphologies, however, can be difficult to interpret. Band- 
shifting can plausibly explain part of the observed changes 
in morphologies; high resolution imaging conducted so far 
mostly in the optical regime traces rest-UV wavelengths at 
Z > 1, which are more affe cted by dust obscuration and on- 
going star formation (e.g., iToft et al.l 120071; lOverzier et al. I 
l2010tlCameron et al. 112010b . Star-forming sites within galax- 
ies may also evolve in size and spatial distribution as a re- 
sult of different dominant triggering mechanisms, and dif- 
ferent physical and dynamical conditions prevailing at ear- 
lier epochs. Kiloparsec-sized "clumps" are in fact ubiqui- 
tous among high-z star-forming galaxies and are believed to 
be at least in part d ue to fragmentation in Toom re-unstable 
gas-rich disks (e.g., Elmegree n et al. II2007L 12009, and refer- 
ences therein). There is now direct kinematic evidence that 
high-z disks are significantly more turbulent (and geometri- 
cally thicker) than their z ~ counterparts, implying accord- 
ingly larger Jeans length setting the characteri stic size scale of 
self-gravitating star-forming complexe s (e.g jBournaud et ail 
120071 120081 iGenzel et al. 112001 l2010h . High gas mass frac- 
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tions at z ~ 1 -2.5 have also been confi rmed through obser- 
vations of CO ro tational transitions (e.g jTacconi et al7ll2010t 
iDaddi et al. 2010). Higher merger rates in the young universe 
may further contribute to the irregular appearance of galax- 
ies. Merger-induced morphological disturbances can how- 
ever last over a wide range of timescales, from < 100 Myr 
up to > 1 Gyr depending on progenitors b aryonic mass ra- 
tio, ga s fraction, and orbital parameters (e.g jLotz et al. 120081 
l2010allbT) . 

Information on the galaxy kinematics is therefore crucial to 
establish the nature and dynamical state of high z galaxies. 
Kinematics provide powerful diagnostics to distinguish reli- 
ably between disks and mergers as a function of the degree 
of sy mmetry of velocity fields and velocity dispersion maps 
(e.g. iKrainovic et al.ll2006l:lFlores et aT1l2006t [Shapiro et al.l 
l2008t lEpinat et al. 2010), and to determine the main source 
of dynamical support from the ratio of rotation velocity to 
local intrinsic velocity dispersion (e.g. iPuech et all 120071 : 
ICresci etal.l 120091: lEpinat et al.l 120091) . Spatially-resolved 
mapping of the two-dimensional kinematics at z > 1 have be- 
come possible during the past decade, with the newest genera- 
tion of sensitive near-IR integral field spectrographs mounted 
on ground-based lOm-class telescopes. 

Using the near-IR integral field spectrometer SINFONI 
at the ESO Very Large Telescope (VLT), we carried out 
the first and largest survey to date of the gas kinematics 
of 80 z ~ 1 — 3 galaxies, "SINS" dForster Schreiber et all 
2006b . 120091: jGenzel et al.ll2006L 120081 : iBouche et al.l 120071: 
Shapi ro et al. 1 12008: Cres ci et aD l2009). The survey includes 
62 rest-UV/optically selected objects at 1.3 < z < 2.6 for 
which we targeted primarily the Ha and [N II] emission 
lines. This "SINS Ha sample" covers roughly two orders 
of magnitude in stellar mass (M+ sa 3 x 10 9 - 3 x 10 11 M ) 
and star formation rate (SFR 10-800 Mgyr" 1 ). The line 
emission is resolved on typical spatial scales of 4-5 kpc 
for the seeing-limited data and down to 1 - 2 kpc for the 
subset of 11 sources followed-up with adaptive optics 
(AO). Kinematic ally, about one-third of the SINS galaxies 
are rotation-dominated yet turbulent disks, another third 
comprises typically more compact dispersion-dominated 
objects, and the remaining sources are interacting/merging 
systems. Concurring results on the kine matic mix at z ~ 1 — 3 
have b een found by other g r oups (e.g..lErb et al.ll2003L|2004j 
12005 IWright etal.l I2007L 120091: itawet al.l l2007bL 120091: 
IStark et al. 1 12001 lEpinat eTaD 120091: Uones et al. I l2010h . If 
the relative proportions vary somewhat between different 
studies, possibly reflecting in part differences in samples 
and mass/redshift ranges, the existence of disks in the early 
universe and their typically high intrinsic local velocity 
dispersion now seem fairly well established. Under these 
conditions, and in the absence of major mergers, internal dy- 
namical processes can drive the secular evolution of disks and 
formation of bulges/spheroids on timescales ~ 1 Gyr, about 
an order of magnitude faster than in pr esent-day disks (e.g . 
Noguchil 119991: llmmeli et al.l l2004allbt ICarollo et all 120071: 
Bournaud. Elmegreen. & E lmegreen 2007; iGenzel et al. I 
2008T lElmegre en et al. 1120091 

While Ha emission provides an excellent probe of the sites 
of ongoing star formation activity and the current dynamical 
state of a galaxy, broad-band rest-frame optical imaging with 
high spatial resolution is required to characterize the struc- 
tural properties of the underlying stellar component. Dynam- 
ical processes at play in the past are imprinted in the stellar 



mass distribution and the stellar population contains a record 
of earlier star formation. Unique insights into early galaxy 
evolution can therefore be gained from the powerful combi- 
nation of ionized gas kinematics and spatial distribution from 
Ha with detailed morphologies from high resolution broad- 
band imaging. 

In this paper, we present the results of deep high resolu- 
tion imaging at 1 .6 fim with the NICMOS/NIC2 camera on- 
board the Hubble Space Telescope (HST) of six of the SINS 
galaxies at 2.1 < z < 2.5. The NICMOS imaging probes the 
rest-frame 5000 A continuum of the galaxies, redwards of the 
Balmer/4000A break. In this range, light from longer-lived, 
lower-mass stars that make up the bulk of the stellar mass be- 
comes more important compared to that of young hot stars 
shining more brightly at shorter wavelengths and the effects 
of dust extinction are smaller than in the rest-UV. Five of 
the targets exhibit compelling signatures of disk-like kine- 
matics. The kinematics of the sixth source unambiguously 
reveal a major merger, although this is not obvious from the 
surface brightness distribution of Ha at seeing-limited reso- 
lution. The original motivation for our NIC2 follow-up was 
to determine whether the stellar light/mass distribution of the 
kinematically-identified disks also follows a disk-like profile, 
and to unveil the components of the merger system. 

With the deep ~ 1 kpc resolution NIC2 1.6 /im imaging, 
our primary goals in the present study are to derive the global 
structural parameters and the detailed morphologies of the 
galaxies as traced by the stellar light. We further explore the 
connection with stellar population properties and, taking ad- 
vantage of the high quality SINFONI data of our targets, we 
test if and how the structure relates to the dynamical proper- 
ties. Our targets extend up to the high-mass end of the rest-UV 
selected z ~ 2 population, and overlap in mass and redshift 
range with other massive galaxy samples with deep and high 
resoluti on NIC2 1 .6 a m data studied b y Ivan )okkum et al.l 
(Hop), iDasyra etal.l d2008h . and iKriek et al.l ( 120091) . The 
similar NIC2 data sets between these samples and ours al- 
low consistent comparisons of the structural and morpho- 
logical properties of massive ?»2- 2.5 galaxies selected 
through different criteria and, anticipating our results, show 
how galaxies of similar masses can be very different in 
those respects, reflecting a wi de range in star formation ac- 
tivity. In a compani on paper (IN. M. Forster Schreiber et aD 
1201 1. in pre paration, hereafter Paper II), we focus on the 
properties of the kpc-scale clumps identified in the NIC2 im- 
ages of our targets, and derive constraints on the nature and 
evolution of these clumps. 

The paper is organized as follows. We introduce the sample 
and the SINFONI data sets, and describe the NIC2 observa- 
tions and data reduction in § [2] In § [3] we revisit the global 
photometric and stellar populations properties of the galaxies 
with the inclusion of the new NIC2 1 .6 /im data, and place 
them in the broader context of z ~ 2 populations in terms of 
stellar mass, star formation rate, and kinematics. We present 
the structural and morphological analysis in §@] and combine 
the NIC2 broad-band imaging together with the emission line 
data from SINFONI in § We compare our sample to other 
massive galaxy samples at z ~ 2-2.6 in terms of structural 
and morphological properties, and address the connection 
with stellar populations in §|6] In §[7] we constrain the spatial 
distribution of stellar mass in one of our targets with addi- 
tional high resolution data from HST/ ACS, compare the rest- 
frame optical and Ha sizes, and discuss the relation between 
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morphologies and Ha kinematics. We summarize the paper in 
§[8] Throughout the paper, we assume a A-dominated cosmol- 
ogy with H = 70h 70 kms" 1 Mpc" 1 , Q m = 0.3, and D, A = 0.7. 
For this cosmology, 1" corresponds to m 8.2 kpc at z = 2.2. 
Magnitudes are given in the AB photometric system unless 
otherwise specified. 

2. SAMPLE, OBSERVATIONS, DATA REDUCTION 

2.1. Choice of NICMOS/NIC2 Targets 

Table Q] lists the galaxies observed along with their photo- 
metric properties. They were drawn from the initial sample 
of 17 rest-UV selected objects at z ~ 2 from the SINS survey 
dForster Schreiber et al.|l2006bl 120091) carried ou t with SLN- 
FONI dEisenhauer et al.l2003HBonnet et al. 12004 at the ESO 
VLT. The galaxies were originally part of the large optical 
spectroscopic survey of z ~ 1.5 -2.5 candidates sel ected by 
their U„GTZ colors desc ribed by lSteidel et al. I (12004 see also 
lAdelberger et al.l 12004 . Additional multi- wavelength data 
include ground-based near-IR J and K s imaging, and space- 
based Spitzer mid-IR photometry at 3 - 8 /im with IRAC and 
at 24 i im with MIPS for th e maj ority of our targets dErb et al.l 
120061 iReddv et al.l 120 101 and IForster Schreiber et al. 1 120091 
for the K s band photometry of MD 41). Prior to the SINFONI 
observations, near-IR long-slit spectroscopy was obtained for 
all of our NIC2 targets with NI RSPEC at the Keck II tele- 
scope and ISAAC at the VLT bv lErbet al.l d200l 120061) and 
Sha plev et al.l d200l . 

The choice of our NIC2 targets was primarily driven by 
their kinematic nature along with their high signal-to-noise 
ratio (S/N), high quality SINFONI data mapping the spa- 
tial distribution and relative gas motions from Ha out to 
radii > 10 kpc (IForster Schreiber et al.l l2006bl) . They in- 
clude five rotating disks (BX663, MD41, BX389, BX610, 
and BX482), and one major merger (BX528). The quantita- 
tive kinematic classificati on performed through application of 
kinemetry is described bv lShapiro et al.l d2008l), and detailed 
dynam ical modeling of the disk s is presented bv lGenzel et al.l 
(2008) and lCrescietal.l d20Q9t). The disks show the expected 
signatures of a fairly ordered two-dimensional velocity field, 
with a monotonic gradient across the source and flattening of 
the velocity curve at large radii, alignment of the Ha kine- 
matic and morphological major axes, and peak of the velocity 
dispersion close to the dynamical and geometrical center. In 
contrast, BX 528 exhibits a reversal in velocity gradient along 
the major axis, as observed in nearby counter-rotating binary 
mergers. Interestingly, the Ha morphology of this merger 
from the seeing-limited SINFONI data does not show features 
that would distinguish it from the disks; all sources have ex- 
tended, more or less irregular/clumpy Ha distributions at a 
typical seeing-limited resolution of ss 0."5. 

Two of our targets (BX663 and BX610) are candidate 
AGN hosts on the basis of their observed mid-IR spectral en- 
ergy distributions fr om Spitzer IRAC and MIPS photometry 
(Redd v~et al.l 1201 61 and N. Reddy, private communication). 
BX 663 also exhibits spectral signatures of Typ e 2 AGN in its 
integrated rest-UV and rest-optical spectrum (Shaple v et al.l 
12004 lErbet al.ll2006h . In our spatially-resolved SINFONI 
data (as well as the new NIC2 imaging presented in this pa- 
per), the host galaxy is well detected and can be distinguished 
from the central compact emission peak with higher [N II]/Ha 
and broad Ha velocity component associated with the AGN, 
underneath the narrower component dominated by star for- 
mation. In contrast, BX610 exhibits no sign of AGN in 



the rest-frame UV and optical, including the spatially- and 
spectrally-resolved SINFONI data, so that its AGN is likely 
to be very obscured and will not affect any aspect of our anal- 
ysis. Therefore, in the context of this paper, we will only 
consider BX 663 explicitly as an AGN source. 

2.2. SINFONI Data Sets 

The SIN FONI data sets of our NIC2 t argets are fully de- 
scribed by Ffars ter Schreiber et al. I d2009l) . They are among 
the deepest of the SINS survey, with on-source integration 
times in the ^T-band (targeting the Ha and [N II] AA6548, 6584 
emission lines) ranging from 3 to 7.3 hours. All galaxies were 
first observed using the largest pixel scale (0."125 pixel" 1 ), 
five of them in seeing-limited mode under very good see- 
ing conditions, and one of them (BX663) with the aid of 
AO using a nearby natural guide star. The typical effec- 
tive angular resolution of these SINFONI /^-band data has a 
FWHM ss 0"5, corresponding to « 4.1 kpc at the redshift of 
the sources (FWHM = 0."39 or 3.2 kpc for the NGS-AO data 
of BX663). Seeing-limited H band data were taken for all 
but BX528 to map the [O III] AA4959,5007 and H/3 emis- 
sion lines, also under very good near-IR seeing of « 0"55 
and using the 0"125 pixel" 1 scale. In addition, Laser Guide 
Star (LGS) AO-assisted K band observations were obtained 
for BX 482 using the intermediate pixel scale of 0." 05 pixel" 1 , 
with total on-source integration time of 6.8 hours and effec- 
tive resolution of 0"17 (1.4 kpc), very similar to that of the 
NIC2 1 .6 /im imaging. 

The central averaged Ha surface brightness (within the Ha 
half-light radius) of the NIC2 targets range from w 1 x 10~ 17 
to w 1 x 10~ 16 ergs" 1 cm " 2 arcsec" 2 . Using the conversion fac- 
tor of lKennicuttl d 1998b between Ha lumino sity and s t ar for- 
mat ion rate, divided by 1 .7 to scale from a ISalpeteri d 19551) 
to a lChabrieU d2003l) IMF, and accounting for beam-smearing 
in computing the intrinsic physical area, we calculate star 
formation rates per unit area of ss 0.15-0.4 MQyr _1 kpc" 2 
(uncorrected for dust extinction). The median S/N per pixel 
of the Ha maps over the regions where line emission is de- 
tected (i.e., S/N > 3) are « 10-15 for the six galaxies. From 
the velocity-integrated line maps and from the noise proper- 
ties o f the data cubes (see details in Forster Schreib er et al.l 
2009), the 3 a Ha surface brightness sensitivities are typically 
« 10~ 17 ergs" 1 cm" 2 arcsec" 2 , corresponding to a limiting star 
formation rate of « 0.03 M Q yr" 1 kpc" 2 . 

2.3. NICMOS/NIC2 Observations 

The NICMOS observations were carried out between 2007 
April and 2007 September with the NIC2 camera onboard 
the HST and using the F160W filter (hereafter referred to as 
the //]6o bandpass). The //i6o bandpass, centered at 1.6 pm 
(FWHM = 0.4 /im), probes the longest wavelengths at which 
the HST thermal emission is unimportant, taking full advan- 
tage of the lack of sky background that limits the sensitivity 
of ground-based near-IR observations. NIC2 has a pixel scale 
of 0"075, critically sampling the HST point spread function 
(PSF) at 1.6 pm, and a field of view of 19."2 x 19"2. While 
the NIC3 camera has greater surface brightness sensitivity, 
NIC2 offers the best combination of resolution and sensitivity 
for our science goals, and also allows more consistent com- 
parisons with the available HST I ACS imaging of MD41 and 
SINFONI+AO observations of BX 482, which have compara- 
ble or better resolution than NIC2. 

Each target was observed for four orbits, with each orbit 
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split into four exposures with a subpixel dither pattern to en- 
sure good sampling of the PSF and minimize the impact of 
hot/cold bad pixels and other such artifacts (e.g., "grot," elec- 
tronic bars). The individual exposure time was 640 s, giv- 
ing a total on-source integration time of 10240 s. The read- 
out mode employed was the SPARS-64 MULTIACCUM se- 
quence with NSAMP =12. We used a square four-point dither 
pattern with step size of 1 "5375 (20.5 pixels). We defined the 
starting positions of each dither sequence as the mid-points of 
the sides of a 0"9 x 0"45 rectangle (12 x 6 pixels). Table 
summarizes the log of the observations. 

2.4. NICMOS/NIC2 Data reduction 

We started with the pipeline-reduced data products, pro- 
cessed with the STSDAS task calnica, version 4.1.1, of the 
memos package within the IRAF environment. We refer to 
the NICMOS Data Handbook^ for details and the descrip- 
tion of specific instrumental features of NICMOS. The steps 
performed in this stage include bias, zero read signal, and dark 
current subtraction, data linearization, correction for so-called 
"electronic bars," flat-fielding, and conversion to count rates. 
The task calnica also identifies cosmic ray hits and various 
bad pixels. These identifications missed on a significant num- 
ber of obvious bad pixels present in our data, so we extended 
the bad pixel lists for each exposure as described below. 

We then used the task pedsky in the same IRAF package 
to subtract a first-order estimate of the background and the 
quadrant-dependent residual bias (or "pedestal") in each indi- 
vidual exposure. We subtracted the median value along each 
column to correct for bias jumps. We constructed a static bad 
pixel mask for each galaxy based on deviant pixels in the nor- 
malized flat field and those identified by calnica with bad data 
quality flags. To identify remaining bad pixels and cosmic 
rays, we first drizzled the images for each galaxy onto a com- 
mon finer pixel grid accounting for their relative shifts, the 
(small) NIC2 distortion, and the static bad pixel masks, us- 
ing the IRAF task drizzle in the STSDAS analysis package. 
A clean reference image was generated by median-combining 
the registered drizzled images with "minmax" rejection. This 
image was inverse-drizzled to the original pixel scale and co- 
ordinate system of each individual exposure with the task blot 
in the STSDAS analysis package. We compared the "blotted" 
reference images with their corresponding input image and 
updated the static bad pixel mask to make a mask for each 
exposure. The background-subtracted exposures were again 
drizzled, now using the updated individual bad pixel masks, 
to create the combined drizzled image of each galaxy. For the 
drizzling, we chose an output pixel scale of 0."05, correspond- 
ing to that of the SINFONI+AO data of BX482, and used a 
square kernel with 0.7 times the size of the input pixels (i.e., 
a "dropsize" of 0.7). We found from experimentation that this 
combination of parameters gave the best results in terms of 
preserving the angular resolution and minimizing the impact 
on the pixel-to-pixel rms noise level of our NIC2 data. 

The combined drizzled images showed substantial residual 
features on scales of a few tens of pixels together with residu- 
als from vignetting at the bottom of the array. These features 
were successfully modeled and subtracted from the combined 
images using background ma ps generated with the S Extrac- 
tor software, version 2.2.2 (Berli n & Arnoutslll996l) . Since 
our science goals require detection of low surface brightness 
emission of the galaxies to large radii, we took particular care 

11 http://www.stsci.edu/hst/nicmos/documents/handbooks/DataHandbookv8 



in making these background maps. We masked out suffi- 
ciently large areas around each object, and then ran SExtrac- 
tor with a background mesh size of 24 pixels and filter size 
of three meshes. This combination provided the best match to 
the spatial frequency of the residual features, without compro- 
mising the faint extended emission from the galaxies. Further 
linear residual features were removed by fitting second-order 
polynomials (with object masking) along columns and rows. 
This step also resulted in a small reduction of the pixel-to- 
pixel rms noise. 

Weight maps (for use in the determination of the noise prop- 
erties and in the structural analysis) were generated in the final 
combination of the exposures with the task drizzle. These re- 
flect the effective total integration time at each output pixel in 
the drizzling procedure (i.e., accounting for the dropsize and 
the relative size of the input and output pixels as well as the 
bad pixel and cosmic ray masking). We normalized the weight 
maps of each galaxy to a maximum of unity. We cropped the 
final images and weight maps to a region 19. "55 x 18 "05 to 
exclude the edges with relative weight < 0.4. 

2.5. PSF 

For the morphological analysis, accurate knowledge of the 
PSF in our data is important. We empirically determined the 
resulting PSF based on stellar profiles in our final reduced 
images. In total, there are four suitable stars (isolated, suffi- 
ciently bright but unsaturated) in three of our NIC2 pointings. 
Their profiles have a Gaussian core and show Airy wings at 
larger radii (around 0"23). To increase the S/N, we averaged 
the peak-normalized stellar profiles. From a Gaussian fit to 
the resulting empirical PSF, the FWHM is 0."145 with ellip- 
ticity of 0.08 (or FWHM of 0."143 and identical ellipticity of 
0.08 for a Moffat fit). Relative to the average PSF, the devi- 
ations of the radial profiles of the stars are at most 15% (on 
average « 5%) out to a radius of 0."35, enclosing 95% of the 
total flux and where the profiles are at w 1 % of the peak value. 
The curves-of-growth of the individual stars differ by < 10% 
(on average w 3%) from that of the average PSF. We adopted 
this empirical average PSF for our analysis. 

2.6. Noise Properties and Limiting Depths 

The quality of the final reduced data is very uniform among 
our targets. To characterize the effective noise properties, 
we fol lowe d the approach described in detai l bv lLabbe et al.l 
(120031) and iForster Schreiber et al. I (|2006al) . For uncorre- 
cted Gaussian noise, the effective standard deviation from the 
mean of the background for an aperture of area A is simply 
the pixel-to-pixel rms er scaled by the linear size N = Va of 
the aperture, a(N) = Na. However, instrumental features, the 
data reduction, and the drizzling and combination procedures 
have added significant systematics and correlated noise in the 
final data. We thus derived the function a(N) directly from 
the data by measuring the fluxes of w 240 non-overlapping 
circular apertures. We placed the apertures at random in 
the maps but so as to avoid pixels associated with sources 
down to S/N sa 3, corresponding to a surface brightness of 
M#i60,ab) = 23.4 mag arcsec" 2 . Series of measurements were 
repeated with different aperture diameters ranging from 0."1 
tol". 

For each aperture size, the distribution of fluxes is symmet- 
ric around the background value (about 0), and is well repro- 
duced with a Gaussian profile. We determined the background 
rms variations from the width of the best-fit Gaussian to the 
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observed distribution for each aperture size. While at any 
fixed spatial scale, the noise properties over the images are 
consistent with being Gaussian, the variations with N deviate 
appreciably from the scaling expected for Gaussian noise, be- 
ing systematically larger and increasing non-linearly. The real 
noise behaviour in our NIC2 images can be approximated by 
the polynomial model: 

NW(a + bN) 
a(N)= (1) 

V w 

where the weight term 1 /y/w accounts for the spatial varia- 
tions in noise level related to the exposure time and is taken 
from the weight maps. The coefficient a represents the ef- 
fects of correlated noise, which are dominated by the resam- 
pling of the images to the final pixel scale through the driz- 
zling procedure. The curvature in the a(N) relation quanti- 
fied by the coefficient b indicates a noise contribution that 
becomes increasingly important on larger scales. For our 
data, this probably originates mostly from instrumental fea- 
tures, and residuals from the background subtraction and flat- 
fielding. Qualitatively, the same second-order polynomial be- 
havior is seen in ground- and space-based optica l and near- 
IR ima ging data analyzed by [Labbe et all (120031 their Fig- 
ure 4). |For ster Schreiber et al. I d200 6a, their Figure 6), and 
IWuvts et al. I (120081 their Figure 3). 

The noise properties are very similar in all NIC2 maps. The 
derived pixel-to-pixel rms values are identical within 3%. The 
mean and standard deviation of the best-fit polynomial coef- 
ficients among the six images are 1.20 and 0.15 for a, and 
0.039 and 0.017 for b. Limiting magnitudes and photomet- 
ric uncertainties from our NIC2 data are computed through 
Eq. [T]for the corresponding aperture size. The effective 3cr 
magnitude limits in a "point-source aperture" with diameter 
d = 1.5 X FWHM = 0."22 (maximizing the S/N of photometric 
measurements in unweighted circular apertures of point-like 
sources) are « 28.1 mag (see Tabled. 

3. GLOBAL PROPERTIES OF THE GALAXIES 
3.1. General Appearance of the NIC2 Images 

Figure [Upresents the final reduced NIC2 images at a sam- 
pling of 0"05 pixel" 1 . The angular resolution from the effec- 
tive PSF with FWHM = 0"145 corresponds to a spatial reso- 
lution of 1 .2 kpc at the median z = 2.2 of the targets. We com- 
puted the geometric center as the unweighted average of the x 
and y coordinates of pixels with S/N > 3. Emission from the 
galaxies is detected out to radii of 0."8 - 1 "2, or 6.5 - 10 kpc. 
The averaged surface brightness of the sources over the emit- 
ting regions (with S/N > 3) is in the range /i(Hi6o.AB) = 
22.6-23.0 mag arcsec" 2 , with typical S/N = 4-6 per pixel. 

The spatial distribution of the Hm,o band light in the images 
is characterized by a diffuse and low surface brightness com- 
ponent over which rich, brighter substructure is superposed. 
The most prominent features are the bright and compact re- 
gions, seen in all galaxies; each of the non-AGN disks con- 
tains several such clumps. BX663 exhibits a bright central 
peak, and only one other and comparatively fainter clump to 
the southeast; since the rest-UV and optical emission lines 
of this object indicate the presence of a Type 2 AGN, the 
compact central H 16 q band peak may not be dominated by 
AGN emission and could trace in large part stellar continuum 
light. Detailed modeling of the optical to mid-IR spectral en- 
ergy distribution of this galaxy will help t o place constraints 
on the contribution of the stars and AGN dK. Hainline et al. I 



1201 1. in pre paration). In BX 528, the merger components are 
clearly separated in our NIC2 image, with a projected dis- 
tance of ss 8 kpc between the compact northwestern and the 
more extended southeastern parts (hereafter BX 528-NW and 
BX528-SE) and significant emission along a "bridge" be- 
tween the two components. The fainter source ss 1"8 to the 
northeast is undetected in our SINFONI data so that its asso- 
ciation with BX 528 is unclear. Likewise for the source ps 2" 
south of BX 482, which could also result from chance align- 
ment. In contrast, the small source pa 0"6 south of the main 
body of BX 389 corresponds to an extension in our seeing- 
limited Ha maps (with about 3 times coarser resolution) and 
is thus likely to be a small satellite at the same redshift and 
projected distance of 5 kpc (hereafter "BX389-S"; see §§|4] 
and|7]). No other source is detected in the NIC2 images within 
ps 2" or projected ps 15 kpc of our primary targets; this area 
corresponds to the deeper part of the SINFONI data, where 
we can search for line emission to faint levels of candidate 
companions to determine a physical association. 

Other characteristics of note include the broad curved fea- 
tures reminiscent of spiral arms seen in the outer isophotes of 
BX 610, on the northwest and southeast sides. In addition, the 
inner isophotes in this galaxy show a significant and system- 
atic increase in ellipticity and change in position angle. From 
fitting elliptical isophotes to the NIC2 image using the task 
ellipse in IRAF, the apparent (uncorrected for beam smear- 
ing) ellipticity and RA. vary from ~ 0.8 and ~ 35° at radii 
r < 0"5 (or 4 kpc) to - 0.1 and ~ 0° at r ~ 1" (8 kpc). The 
interpretation of the inner isophotal feature is complicated by 
the presence of the clumps, which are mostly aligned along a 
similar R A., but it is suggestive of a bar-like structure as seen 
in local barred spirals. 

3.2. Photometric Properties 

The new Hia) band photometry from our //S7YNICMOS 
data constitute a significant addition to that available for our 
targets, especially for BX482 without any previous near-IR 
imaging and MD41 for which the existing /f-band data are 
fairly shallow. Published data include ground-based optical 
U„G1Z imaging for al l sources, and ne ar-IR JK * imaging for 
a subset, presented by Erb et al. ( 2006, see also Steidel et al. 
2004). The photometry was measured in matched isophotal 
apertures determined from the 1Z band images, and corrected 
to "total" magnitudes in elliptical apertures scal ed based on 
the fi rst moment of the 7^-band light profiles ( Steid eTet al. I 
2003). Spitzer/IRAC mid-IR photometry has subsequently 
been obtained for all objects except BX482 (Redd y et al. I 
2010). To construct the broad-band spectral energy distribu- 
tions (SEDs) and derive the global stellar and extinction prop- 
erties of our targets, we however excluded the IRAC data for 
the following reasons. The optical and near-IR photometry 
provides the most uniform data sets among the six galaxies 
in terms of wavelength coverage whereas of the five objects 
with IRAC data, only three were observed in all four channels. 
Moreover, contamination by AGN emission in BX 663 and 
BX 610 becomes important at the longer wavelengths probed 
by IRAC and would affect the modeling results. 

To include the Hm,o band in the SEDs of our targets, we 
measured the total magnitude in a circular aperture of diam- 
eter 3" on the NIC2 images (indicated in FigureQ]). Because 
of the small NIC2 field of view, the accuracy of the measure- 
ments in the 7^-band isophotal aperture used for the ground- 
based photometry on the PSF-matched //i6o band images is 
compromised by the difficulty of determining the local back- 
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ground sufficiently reliably: the area empty of sources signal 
after convolving to the optical seeing of ~ 1" is too small. 
The median difference between the total magnitudes based on 
the fixed circular aperture measurement on unconvolved NIC2 
maps and that in the isophotal apertures on PSF-matched 
NIC2 maps is -0.16 mag, with the maximum for BX389 
(-0.51 mag) and best agreement for BX528 (-0.002 mag). 
The photometry is given in Table Q] and the optical/near-IR 
SEDs are plotted in Figure|2] 

By their primary selection, all targets s atisfy the U„G1Z 
color criteria of "BX" o r "MD" objects (Adelberge r et al.l 
I2004t iSteidel et al. 1120041) . Synthetic broad-band colors mea- 
sured on the spectrum of the best-fitting st ellar population 
model to the observed SEDs (described in § !3.3l below) sug- 
gest all targets would also satisfy the BzK = (z-K)p&-(B- 
z)ab > -0.2 mag crit erion of z> L4 "s tar-forming BzK ob- 
jects" introduced by Daddi et al. (2004). In addition, two 
of the four sources with both J and K photometry (BX 663 
and BX 389) have the r ed (J- K S ) AB > 1 .34 colors of "Distant 
Red G alaxies" (DRGs; iFranx et al. 1120031; Ivan Dokkum et al. I 
120041) . and one (BX 610, (J-K S ) AB = 1 -28 mag) is within 0.3 a 
of the DRG criterion. For these three sources and for BX 528, 
the photometry is consistent with the presence of a significant 
population of evolved stars along with the young stars produc- 
ing the blue rest-UV SED. The best-fit models to the SEDs of 
these four objects imply rest-frame (U-B)ab ~ 0.75 -0.9 mag 
colors, clos e to or just at the separation adopted by, e.g., 
iRriek et aTTI (12009b . between "red" and "blue" objects. In 
contrast, MD41 and BX482 are significantly bluer with rest- 
frame (U-B)ab ~ 0.5 mag. 

Our NIC2 sample is relatively bright in the near-IR; the 
five targets with K band imaging have a median K s ^b = 
21 .8 mag and span a range of w 1 .2 mag in observed K band 
(^s,ab = 21.1-22.3 mag) as well as in rest-frame absolute V 
band magnitude (as derived from the SED modeling below, 
with Mv.ab between -23.1 and -21.9 mag). In addition, the 
Spitzer/MlPS 24 /im measurements available for four of the 
targets (BX389, BX610, BX528, and BX663) indicate that 
they are among the brightest 3% of rest-UV s elected z ~ 2 
galax ies, in terms of rest- frame 8 /im luminosity (Redd vet al.l 
l2O"l0l) . 

3.3. SED Modeling and Stellar Properties 

Stellar population synthesis modeling of BX528, BX663 , 
BX389, and BX 610 was first presented bv lErb et al.l (I2006T) . 
The new H\(,q photometry from NIC2, and K band flux 
for MD41, provide additional constraints and allow us to 
extend the SED modeling to BX482 and MD41, which 
had previously only three optical data points. For consis- 
tently derived stellar and dust extinction properties within 
our sample, we re-modeled the optical/near-IR SEDs of 
all six objects^ We followed standard procedures (e.g. 
Shapley_etaL][2001; Erb et al.ll2006t IForster Schreiber et al. I 
2006al) . exploring also a range of parameters and model in- 
gredients to assess possible systematics from the assumptions 
adopted. Our specific proce dure is described in detail by 
Forster Schreiber et al. (2009) and summarized here. 

We used the stel lar evolutionary code from 
iBruzual & Charlotl ( 120031) . The age, extinction, and lu- 
minosity scaling of the synthetic model SED were the 

12 We verified that the SED modeling results are essentially unchanged 
when including the available IRAC photometry for the three targets without 
AGN contribution to their mid-IR SED. 



free parameters in the fits. The ages were allowed to vary 
between 50 Myr and the age of the Universe at the redshift 
of the source. The lower limit is set to avoid implausibly 
young solutions and corresponds approximately to the 
dynamical timescale Td yn = rj v from the inclination-corrected 
velocity v at a radius of r ~ 1 kpc of the galaxi es (e.g. 
IForster Schreiber et al.l l2006bl: iGenzel et all 120081) . We 
e mployed the | Chabrierl ((2003) IMF and the reddening curve 
of Calzetti et al. (2000), and assumed a solar metallicity. We 
fixed the star formation history (SFH), adopting a constant 
star formation rate ("CSF") model, which we found appro - 
priate for most galaxies (see below, and also Erb et al. 2006). 
We derived the formal (random) uncertainties on the best-fit 
parameters from 200 Monte-Carlo simulations, where we 
varied the input photometry assuming a Gaussian proba bility 
distribution for the measurements uncertainties (see § 12.61 ). 
To explore the impact of model assumptions and ingredients, 
we als o used the evolutionary synthesis code of M arastonl 
(120051) . and ran suites of models with exponentially declining 
star formation rates (SFRs) with e-folding timescales r 
between 10 Myr and 1 Gyr, with metallicity of 1/5 and 2.5 
times solar, and with ex tinction laws appropriate for the 
Small Magellanic Cloud (iPrevot et al.l 1 1984 iBouchet et"aT1 
1985) and the Milky Way dAllenlll976l) . 

Table [3] lists the best-f it properties for our adopted set of 
IBruzual & Chariot (2003) models with CSF, solar metallicity, 
and lCalzetti et al.l (12000) law along with the formal uncertain- 
ties taken as the 68% confidence intervals from the Monte- 
Carlo simulations around the best-fit values. Figure [2] shows 
the co rresponding synthetic model spectra; best-fit Maraston 
(2005) models for the same CSF, metallicity, and reddening 
law are also shown for comparison. To gauge systematic 
uncertainties, we looked at the variations in derived proper- 
ties for the other suites of model assumptions and synthesis 
code considered. For individual galaxies, and compared to 
the adopted properties, the median of the variations are as fol- 
lows: the ages are mostly younger by a factor of ss 2, the 
extinction in V band Ay decrease by « 0.2 mag, the stellar 
masses tend to be lower by a factor of w 1 .5, the absolute 
and specific SFRs tend to decrease by factors of « 2. The 
overall sense of these variations reflects in part the fact that, 
by construction, CSF models generally lead to older ages than 
declining SFHs. While the systematic uncertainties above are 
significant, the derived properties of the galaxies are robust in 
a relative sense. In particular, MD41 and BX482 remain the 
youngest and least massive of the sample in any case. 

Our choice of CSF is obviously very simplistic. For the 
range of SFHs above, we find that the longest timescales (r = 
1 Gyr or CSF) are favored (in a minimum \ 2 sense) for four 
out of six galaxies. MD41 and BX482, with bluest SEDs, 
have best fits formally obtained for models with r < 50 Myr 
but the age, M + , and SFR are within 3 a or less of the values 
for CSF. Recent work has emphasized the fact that declining 
SF Hs may well be inappropriate for high z star-forming galax- 
ies ([Renzini 2009; Mara ston et al.ll20lol) . Our exploration 
of constant and declining SFHs was motivated by continuity 
with previous work and consistency with the SED modeling 
assumptions for the comparison samples discussed in § [6j 

We did not account for emission line contribution in our 
SED modeling. The main contribution expected for our sam- 
ple is from Ha in the K band, and these are in the range 
7% -19% (with median of 1 3%) based on o ur SINFONI K 
band observations (see also lErb et al.l 120061) . The median 
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[N II] A6584/Ha ratio is 0.23 (with highest ratio 0.43 for 
BX 610). SINFONI H band data are available for five targets; 
the [O III] doublet and H/3 are all detected in two cases and 
their contributions to the Hk,q flux densities are smaller than 
for Ha in K, and 3 a upper limits for the others are at most 
29% (see § |5j. The Ha-corrected K band flux densities are 
indicated in Figure|2l the impact on the SED modeling results 
is not significant in view of other uncertainties involved. 

Our targets span over an order of magnitude in stellar mass, 
from 7.7 x 10 9 to 1.0 x 10 n M Q with median 5.3 x 10 10 M©. 
The stellar populations that dominate their rest-UV to opti- 
cal SEDs cover the full range allowed in our modeling (from 
50 Myr to 2.75 Gyr), and are moderately obscured (with 
Ay between 0.6 and 1.2 mag, or equivalently E(B-V) = 
0.15-0.30 mag). Four of the six targets have in fact ages 
consistent with the age of the Universe at their redshift, indi- 
cating the presence of mature stellar populations. The median 
absolute and specific SFRs are m 50 M©yr _1 and 0.6 Gyr" 1 
(with ranges of 25-185 M Q yr _I and 0.6-24 Gyr -1 ). Fig- 
ure^ compares our NIC2 sample in the M t - SFR plane with 
the full SINS Ha sample, and a purely /f-selected sample 
taken from the Chandra Deep Field South (CDFS) in the same 
z= 1.3-2.6 interval and to the same A" s v apa = 22 mag limit 
as the SINS sample ( For ster Schrei ber et al. 2009, and also 
IWuyts et al.ll2008l) . The properties of these comparison sam- 
ples were derived using the same SED fitting procedure and 
model assumptions as for the NIC2 sample. Clearly, and by 
selection, our NIC2 galaxies probe the actively star-forming 
part of the z~2 galaxy population (see also discussions by 
lErbet al.l l2006: Fors ter Schreiber et al.ll2009l) . Their median 
stellar mass and SFR are comparable to those of the SINS Ha 
sample («3x 10 10 M© and 70 Moyr" 1 ). 

3.4. Ha and Kinematic Properties 

If our targets do not particularly stand out in terms of their 
stellar properties compared to the SINS Ha sample, they 
do more so in terms of their Ha brightness and kinematics. 
This distinction arises because, as described in §[2] we chose 
them among the sources with highest S/N and spatially best- 
resolved Ha emission from our seeing-limited SINFONI data 
sets. Five of them were explicitly selected because of their 
disk-like kinematics, and BX 528 because it shows in contrast 
kinematic signatures expected for counter-rotating mergers. 
The spatial extent and brightness of our targets make them 
particularly well suited for a detailed study of morphologies 
and reliable assessment of structural parameters from sensi- 
tive kpc-resolution continuum imaging. 

Figure [3j> compares the NIC2 galaxies with the SINS 
Ha sample in the velocity-size plane, with the half- 
light radii and (in clination-corrected) c i rcular velocit i es de- 
rived from Ha (IGenzel etal. I 12001 ICresci et al. I 120091 
IForster Schreiber et al. 1 120091) . The NIC2 targets lie at the 
high end of the Vd — rirt(Ha) distribution. With ^^(Ha) f=s 
4-5 kpc, they are all larger than the median of 3. 1 kpc for the 
SINS Ha sample. The va for the disks are roughly equal to or 
larger than the median for SINS galaxies of ss 180 kms" 1 , 
and up to ~ 300 kms" 1 . For the merger BX528, the ve- 
locity gradient is lower, and leads to a lower equivalent Vd 
of 145 kms" 1 (see Forster Schre iber et al. 1 120091) . The ra- 
tios of circular velocity to intrinsic local velocity dispersion 
of the disks are in the range Vd/ero ~ 2-6, which is sig- 
nificantly lower than for local spiral galaxies (values in the 
range - 10-40; e.g jDib. Bell. & Burkertl2006l:IEpinat et al.l 



2010) and suggests comparatively larger gas turbulence and 
geometrical thickness. These low Vd/co ratios appear to 
be a char acteristic feature of early disk galaxies a t z 
1 - 3 (e. g. IForster Schreiber et al l|2006bl: IWright et al 1 1 2007 
I Genzel et al.l 120081 IStark et al.l 12008 ICresci et alTl2009 



lEpinatetal.1120091) . 

The median integrated Ha flux of the NIC2 targets (uncor- 
rected for extinction) is F(Ha) = 1 .9 x 10~ 16 erg s" 1 cm" 2 (and 
range from 1.1 to 3.1 x 10~ 16 ergs _1 cm" 2 ). This places them 
among the brighter half of the SINS Ha sample (with median 
F(Ha) = 1.1 x 10~ 16 ergs" 1 cm" 2 ). In terms of integrated Ha 
line widths, their median is er(Ha) = 155 kms" 1 (spanning the 
range 120-245 kms" 1 ), somewhat broader than the median 
for the SINS Ha sample (130 kms" 1 ). This difference is due 
in part to the large observed velocity gradients for the sources. 

4. ANALYSIS OF THE REST-FRAME OPTICAL MORPHOLOGIES 

4.1. Parametric Analysis of the Galaxies' Structure 

4.1.1. Methodology 

In a first step to characterize quantitatively the morpholo- 
gies of our galaxies, we followed a parametric approach. The 
main goal here is to derive the global str uctural parameters of 
the sources. We used the code GALFIT dPeng et al.ll2002l) to 
fit the t wo-d imensional surface brightness distributions with a 
Islrsicl (fl968l) profile: 



/(r) = 7(0) exp 



-b n (r/R e ) 



l/n 



(2) 



where I is the intensity as a function of radius r, b„ is a 
normalization constant, and R e is the effective radius enclos- 
ing half of the total light of the model. In this parametriza- 
tion, ellipticals with de Vaucouleurs profiles have a Sersic in- 
dex n = 4, exponential disks have n = 1, and Gaussians have 
n = 0.5. In the fitting procedure, GALFIT convolves the model 
light distribution to the effective resolution of the data, and de- 
termines the best fit based on x 2 minimization. We input our 
empirically determined PSF, and the weight maps scaled by 
the pixel-to-pixel rms measured on the reduced images (see 
§|2]). To reduce the impact of large-sc ale r esiduals from flat- 
fielding and background subtraction (§ !2.61 l. we performed the 
fits within a region 7" x 7" centered on the sources. 

The free parameters in the fits were the Sersic index (n), 
the effective radius (R e ), the axis ratio (b/a, semi-minor over 
semi-major axis radii), the position angle (RA.) of the ma- 
jor axis, and the total magnitude. Given that surface bright- 
ness fitting of faint galaxies is sensitive to the background 
determination and that Sersic profile para meters tend to be 
degenerate with the "sky " value (see, e.g. JPeng et al. 1120021 ; 
van Dokku m et al. 1 12008). we fixed the sky level S in our fits. 
Moreover, asymmetric and clumpy surface brightness distri- 
butions can strongly bias the results, so we also fixed the coor- 
dinates xq and yo of the center position. In order to derive re- 
alistic uncertainties, we performed a series of fits where these 
fixed parameters were drawn at random from plausible distri- 
butions. Specifically, for every galaxy we ran GALFIT 500 
times, each time varying the xq and yo according to a uni- 
form distribution within a squared bo x of 2x2 pixels around 
the geometric center as defined in § 13.11 which corresponds 
closely to the dynamical center of our targets. For the back- 
ground, the fixed input S values were drawn at random from 
the central 10% of the distribution of individual pixel values 
in the regions empty of sources (with this range reflecting the 
uncertainty in the background level). The best-fit parameters 
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were taken as the median of the results for all 500 runs, and 
the 1 a uncertainties as the 68% confidence intervals about the 
median. 

We performed additional series of fits to explore systemat- 
ics from the center and sky background determinations. In all 
cases, the initial guesses for the center and sky values were 
generated as described above. In the first two series, one of 
the sky or center was chosen as a free parameter and the other 
was fixed. In a third series, both parameters were free. For 
BX663, MD41, BX389, and BX610, the best-fit structural 
parameters for these additional series are all within the 68% 
confidence intervals of the results for the first series with sky 
and center fixed to the randomly drawn values in each iter- 
ation. For BX 528 and BX 482, the two series with the free 
center position lead to significantly different results as ex- 
pected from their particularly asymmetric morphologies, with 
changes in R e and n by > 2 a. 

For the interacting system BX528, single component fits 
are clearly inappropriate. When the center is free to vary, the 
fits are heavily influenced by the brighter BX 528-SE com- 
ponent: the best-fit center moves towards this component, the 
Sersic index increases to match the steeper inner profile and 
the R e also increases. For this system, we thus also ran a 
series fitting two Sersic profiles simultaneously, fixing the re- 
spective center position and sky to random values drawn from 
appropriate distributions as explained above. Likewise, we 
performed two-component fits for BX 389 to derive the struc- 
tural parameters of the small southern companion BX 389-S. 
The best-fit parameters for the main part of BX 389 are al- 
most identical to those obtained with the single-component 
fits. For BX 482, the large and bright clump on the southeast 
side (which makes up w 15% of the total light after subtract- 
ing the "background" light from the host galaxy; see Paper II) 
drives the center in this direction when allowed as a free pa- 
rameter, and leads to a somewhat higher best-fit n index. At 
the same time, because the northwest side of the galaxy is 
overall fainter and lacks bright small-scale features, the best- 
fit R e is smaller. We did not consider two-component fits for 
BX482 because the kinematics as well as morphologies are 
more consistent with a single system. 

4. 1 .2. Results for the NIC2 Targets 

Table|4]gives the results for the main structural parameters, 
R e , n, b/a, and P. A., of our targets from the single-component 
fits as well as from the two-component fits for BX 528 and 
BX389. In Figure|4] the first two columns show the best-fit 
values for single-component fits along with the distributions 
obtained from the 500 GALFIT runs in each parameter. In the 
third column, ellipses are drawn on the NIC2 images showing 
the effective radius, axis ratio, and position angle of the best- 
fit models. The fourth column shows the fit residuals, reveal- 
ing the clumpy structure and large scale asymmetries in the 
light distributions. Figure [5] shows the same but for the two- 
component fits to BX528 and BX389. Figure [6] compares 
the radial light profiles of the galaxies' data and of the best- 
fit model (here using the two-component model for BX 528). 
The profiles correspond to the average along ellipses centered 
on the best-fit xq and yo position, and with best-fit axis ra- 
tio and position angle of the single-component models. This 
tends to smooth out the small-scale structure and asymmetric 
features of the surface brightness distributions; the detailed 
substructure will be considered in the next subsection. 

All of our targets qualify as morphologically late-type, 
disk-dominated systems, often defined as having n < 2-2.5 



fe.g.. iBell et aDl2004t iTruiillo et al. Il2006h . The only excep- 
tions are BX528-NW and BX 389-S that have n w 3 but the 
large uncertainties make them within 1 a of the above cri- 
terion. A remarkable feature of the four non-AGN galaxies 
with disk-like Ha kinematics is their shallow inner light pro- 
files. With Sersic indices in the range « 0.15-0.6, these 
are significantly less centrally concentrated than exponen- 
tial disks (by definition n = 1). The residuals are overall 
smallest for MD41 and BX389, the two most edge-on disks 
based on their axis ratios (b/a « 0.30) from the fits to the 
i/i60 band ima ges as well as from their Ha morphologies 
and kinematics |Fqrster Schreiber et al.H2006bl:lGenzel et al. I 
I2008t ICresci et al. II2009D . The spiral-like features in BX610 
noted in § 13.11 are most obvious in the outer parts beyond 
the effective radius. The residuals for the single- and double- 
component fits to BX 528 indicate that neither case provides 
a very good representation. There is clearly significant light 
emitted along the "bridge" between the SE and NW compo- 
nents. 

BX 663 has the highest n index among our sample, reflect- 
ing the fairly prominent central peak and shallower profile 
of the extended emission at larger radii. We explored two- 
component models, meant to represent the case of a disk + 
bulge and of a disk + unresolved nuclear source. A free n 
component was considered for the disk, and the second com- 
ponent was taken as either a de Vaucouleurs profile with fixed 
n = 4 or as the empirically determined PSF Because the cen- 
tral peak is accounted for by the n = 4 or PSF component, the 
underlying disk ends up with low n of 0.24 (disk + bulge) or 
0.63 (disk + unresolved source), in the range obtained for the 
other five SINS targets. The disk R e values, however, differ 
by 10% or less (< 2.5 a) from the best-fit single-component 
R e . For the disk + bulge case, the de Vaucouleurs component 
has R e = 2.47±g "5 kpc. In terms of \ 2 values and of residu- 
als, the best-fit two-component models are indistinguishable 
from the best-fit single Sersic model, so we consider only the 
single-component model for BX 663 throughout this paper. 

We find no significant correlation between structural pa- 
rameters (R e , n, and b/a) and SED derived properties (M*, 
age, Ay, absolute and specific SFR) among our targets. Stud- 
ies of larger samples based on ground-based seeing-limited 
or higher resolution NIC MOS/NIC 3 near-IR imaging (e.g., 
iTrujillo et al.ll2lMlTbfte't al.ll2007l:lFranx et al. 2008) show 
the existence of relationships between size, stellar mass, and 
star formation activity at z ~ 2 - 3 (though with significant 
scatter). Among our five disks, MD 41 and BX 389 with high- 
est Ay values have the lowest axis ratios, which seems con- 
sistent with the trend of higher extincti on in more edge-on 
star-forming di sks observed locally (e.g.. lMaller et al. Il2009h 
lYip et al.1120 10). However, clearly, only the strongest correla- 
tions would be discerned with our very small sample. 

4.2. Non-parametric Analysis of the Morphologies 
4.2.1. Methodology 

The models considered in the previous section are admit- 
tedly very simplistic. Simple one- or two-component axisym- 
metric models can reproduce reasonably well the overall sur- 
face brightness distribution of some of our sources but do 
poorly for others. Such a parametric approach appears to be 
satisfactory to estimate global parameters such as the effec- 
tive radius and shape of the radial light profile. On the other 
hand, it obviously does not capture the prominent small-scale 
and irregular structure of clumpy disks and merging systems. 
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The difficulty of describing star-forming galaxies at z ~ 2 - 3 
within the framework of a Hubble sequence of more regu- 
lar spiral and elliptical morphological types has been high- 
lighted previously by several author s (e.g.. lLotz et al. 1 120061 : 
lLaw et al.ll2007at IPeter et all 120071) . Accordingly, we also 
followed a non-parametric approach to characterize in detail 
morphologies of our galaxies, measuring the Gini, M20, and 
Multiplicit y coefficients , which have been introduced a nd ap- 
plied by AbrahametaD d2003l) . lLotz et all d2004l2006h . and 
lLaw et all d2007al) . 

The Gini parameter, G, is a measure of the uniformity with 
which the flux of a galaxy is distributed among its constituent 
pixels. G = corresponds to the case where the flux is dis- 
tributed evenly among all pixels, while G = 1 indicates that 
a single pixel contains all of the flux. The Gini parameter is 
formally defined such that: 



G = 



1 



fN pix (N pix -l) 



(3) 



where / is the mean sky-subtracted flux per pixel, /]■ is the 
sky-subtracted flux in pixel 2, and the pixels have been sorted 
in order of increasing flux. 

The M20 parameter represents the second-order moment of 
the brightest 20% of a galaxy's flux. The total second-order 
moment is defined such that: 



Ay, 

i 



-x.y + iyi-y,) 2 ] 



(4) 



where /,■ is the flux in each pixel, x, and indicate the pixel 
location, and x c and y c indicate the centroid of the galaxy light 
distribution, as determined from its first-order moments. M20 
is then computed by sorting galaxy pixels by flux, and sum- 
ming M, until the pixel is reached at which 20% of the total 
galaxy flux is contained in pixels that have brighter or equal 
fluxes. That is: 

M 20 = log 10 , while Y^ft < 0.2f tot . (5) 



lLotz et all (120041) introduced M20 as an indicator sensitive to 
merger signatures such as multiple nuclei. 

Finally, a third c lassificatio n param eter, the Multiplicity 
(<]/), introduced by lLaw et aL (2007a), is somewhat analo- 
gous to the M20 parameter, in that it is sensitive to the presence 
and distribution of multiple clumps of flux. "J is conceptually 
akin to the total gravitational potential energy of a massive 
body. The "potential energy" of the pixel flux distribution can 
be described as 



ipn 



(6) 



where rn is the distance between the zth and jth pixels. The 
greatest "potential energy" for a given set of pixel fluxes is 
the most compact distribution, with the brightest pixels at the 
center, and flux decreasing outwards. Such a configuration 
would be described by 



AV> AW ,, 
JJJ± 
/ j r l ' 

i=l j=lJ¥> ij 



(7) 



where r-y is the distance between the jth and jth pixels in 
the maximally compact configuration. The multiplicity, 'J, 
is then defined to be 



* = lOOlog 
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compact 



(8) 



lLaw et all (l2007al) argued that provides larger dynamic 
range than M2o- Illustrative examples of the relationship be- 
twee n multiplicity and ov erall appearance ar e provided by 
both lPeter et ail (120071) and lLawet all (I2007a[) . 

Some care must be taken in assigning pixels to a galaxy 
light distribution, especially when considering galaxies over 
a range of redshifts, where the differential effects of cosmo- 
logical surface brightness dimming must be taken into ac- 
count. We search for pixels associated with a galaxy within 
a l."5 radius of the galaxy centroid. Pixel selection is per- 
formed on a lightly smoothed version of the data, which has 
been convolved with a circularly symmetric Gaussian kernel 
with a 1 pixel standard deviation. We select pixels to be as- 
sociated with the galaxy if the smoothed flux is more than 

1.5 y ^ x a above the sky, with a calculated from the 

raw data, and where Zhigh is the highest redshift for the subset 
of galaxies under consideration. The redshift-dependent fac- 
tor counteracts the effects of surface brightness dimming in 
a given passband. Once pixels are assigned to a galaxy light 
distribution, the raw, unsmoothed data are used for morpho- 
logical analysis. 

4.2.2. Results for the NIC2 Targets 

Figure [7] shows the resulting Gini, M20, and ^> coefficients 
for the SINFONI sample. To control for surface brightness 
dimming effects, each galaxy was analyzed as if it would be 
observed at the maximum redshift in the sample, i.e. z = 2.43. 
The Gini values range from G = 0.25-0.36, with a median 
of G me d = 0.33. The \P values are in the range 2.8-9.8, with 
a median of ^med = 5.4. Finally, the M20 values range from 
—1.1 to -0.6 with a median of M20.med = -0.7. 

In order to interpret these numbers, we compare the NIC2 
morphologies for objects in our sample with those derived 
from much larger samples of rest-UV selected galaxies at sim- 
ilar redshifts with morphologies from optical imaging with the 
HST Advanced Camera for Surveys (ACS ). The ACS F 814W 
O'814-band) rest-UV morphologies from IPeter et ail d2007l) 
provides one such comparison sample, spanning a compa- 
rable redshift range to that of the objects in our NIC2 sam- 
ple. To make a controlled compari son, we re strict the analy- 
sis to the subset of objects in the IP e t er et al. I sample with the 
same average pixe l -based S/N ratioQ This subset of objects 
in the IPeter et ail (J2007) sample has an identical median G 
value to the objects in our sample, and median ^ and M20 
values that are very similar (^med = 7.0 and M2o,med = -0.8). 
The dataset used to analyze the ACS iR^-band mo rphologies 
of rest-UV selected galaxies by IPeter et al. I (2007) ha s simi- 
lar dep th to the summed BViz images featured in lLaw et al.l 
(2007a), where rest-UV selected objects were also analyzed. 
Accordingly, these authors found roughly identical rest-UV 
morphological results. Overall, the rest-frame optical mor- 
phologies for the objects in our sample are similar to the typi- 

3 To determine the average pixel-based S/N, we calculated the average 
ratio over all object pixels of the flux in each pixel and the pixel-to-pixel noise 
in the background. 



10 



Forster Schreiber et al. 



cal rest-frame UV morphologies for z ~ 2 star-forming galax- 
ies. 

One of the objects in our sample, MD 41, has existing deep 
ACS F814W imaging. In Figure |8] both ACS and NIC2 im- 
ages are shown side-by-side at their respective original spa- 
tial resolution, along with an RGB color composite after PSF 
matching of the A CS m ap to the resolution of the NIC2 map 
(as described in § I7.lt . The average pixel S/N of MD41 is 
very similar in the two images, making for a fair comparison 
of the wavelength-dependent morphology. Overall, the mor- 
phologies in the two bands are strikingly similar. The most 
notable difference consists of the two brightest clumps in the 
ACS map at the southwestern edge of MD41, which have 
only faint counterparts in the NIC2 image and, accordingly, 
exhibit bluer colors than the rest of the galaxy. In addition, 
two compact (radius ~ 0"2) sources offset by ~ 1" from the 
main body of the galaxy appear in the ACS image, alone, 
indicating significantly bluer colors than the main region of 
galaxy emission. We have not yet determined if these sources 
are associated with MD41. In particular, a search of the SIN- 
FONI data cube does not reveal any emission line at these 
locations and within ±1000 tans" 1 of the Ha line of MD41, 
which could support a physical association. 

The Gini values measured in the ACS and NIC2 images 
of MD 41 are both G ~ 0.3, while the multiplicity is slightly 
lower for the rest-frame optical image (6.3) than the rest- 
frame UV image (9.7 with the two offset objects included 
and 8.1 without). The Moq values are comparable as well, 
with -0.7 in the rest-frame optical and -0.9 in the rest-frame 
UV. While the parameters and appearance in each band are 
not identical, we stress that the differences are not signif- 
icant. MD41 does not display a strong morphological "k- 
correction," according to which the rest-frame UV and opti- 
cal morphologies would trace light from very different stellar 
populations. Similar comparative studies of the rest-UV and 
optical morphologies of samples of z ~ 1 .5 -4 galaxies based 
on ACS op tical and recent WFC 3 H k,q band imaging are pr e- 
sented bv lOverzier et al. I d2010l) and lCameron et alT d2010). 

For another object in the sample, BX482, the SINFONI 
Ha map obtained with LGS-AO and 0"17 PS F FWHM 
dGenzel et all 120081: IForster Schreiber et al.ll2009h is of suf- 
ficient resolution and S/N to measure the detailed morpho- 
logical parameters using the same methodology and compare 
with the NIC2 Hi^o band properties. We plot the morpholog- 
ical parameters derived from the Ha map of BX 482 in Fig- 
ure [7] along with the NIC2 -based parameters. The Ha line 
and //i6o band morphologies are quite similar in terms of G, 
^, and M20, with the prominent high surface brightness clump 
appearing in both SINFONI and NIC2 images. 

We searched for correlations between the morphological 
(G, ^P, and M20) and stellar population parameters (M+, age, 
Ay, absolute and specific SFR). The most significant (3 a) 
trend is found between stellar mass and G, such that more 
massive galaxies are characterized by higher G values. Cor- 
relations between G and other gauges of galaxy evolutionary 
state such as age and specific SFR, are not readily apparent 
within our sample. However, our sample of six objects is ob- 
viously too small for discerning all but the strongest trends. 
Furthermore, with only actively star-forming galaxies in our 
sample, we do not probe the full dynamic range of z ~ 2 star- 
forming histories. In section |6j we return to the connection 
between morphology and other galaxy properties, consider- 
ing our SINS sample alongside z ~ 2 galaxies selected using 



very different criteria (|van Dok kum et al. 2008; Dasvr aet al.l 
l200aiKriek et al.ll2009h . 

5. COMPARISON WITH EMISSION LINE PROPERTIES FROM 
SINFONI 

In this section, we combine the NIC2 images and SINFONI 
emission line data of our six targets to compare the Ha mor- 
phologies and structural parameters with those from the rest- 
frame optical emission, and determine the contribution from 
[OIII] and Hf3 line emission to broad-band //160 flux densities. 

5.1. Comparison of Morphologies 

For proper comparison of the H^q band and Ha maps, 
we first convolved the NIC2 images to the spatial resolution 
of the SINFONI data. The SINFONI FOV is too small to 
include neighboring stars observed simultaneously with the 
galaxies, and the PSFs are based on the average of broad- 
band images synthesized from the SINFONI data cubes of 
stars observed immediately before or after the target. These 
non-simultaneous PSF calibrations are however estimated to 
provide good representations o f the PSF of the targets' data, t o 
« 20% in terms of the FWHM ( IForster Schreiber et al.l2"009h . 
We computed the smoothing kernels with a Lucy-Richardson 
deconvolution algorithm. Curve-of-growth analysis of the 
convolved NIC2 and original SINFONI PSFs indicates that 
the fraction of enclosed flux agrees to better than 5% in cir- 
cular apertures with diameters in the range 0"15- l."5 in all 
cases. 

After convolution, we resampled the NIC2 images to 
the pixel scale of the SINFONI data where appropriate 
(Of 125 pixel" 1 for all sets but the AO 0"05 pixel" 1 data of 
BX482) and aligned them with the emission line maps re- 
lying primarily on the outer isophotes. The detailed emission 
line and broad-band morphologies can plausibly differ for any 
galaxy. However, for our sources, and especially at a seeing- 
limited resolution of w 0."5, the overall light distribution turns 
out to be very similar. We verified the registration with line- 
free continuum maps extracted from the SINFONI data cubes 
for the brighter sources. In contrast to the line maps, the con- 
tinuum maps have lower S/N so that mostly only the central 
brighter regions or the continuum peak are reliable. The align- 
ment based on line and continuum in those cases agrees well 
within 1 - 2 SINFONI pixels. 

Figures [9] and [10] show, for each target, the original Hk,q 
band image along with the Ha line map and velocity field 
with contours of the convolved Z/160 band image overplotted. 
The most detailed comparison is possible for BX482 with 
AO-assisted SINFONI data. Both the lower surface bright- 
ness component and the brighter sub-structure are compara- 
ble between the H^q band and Ha emission. This similarity 
is reflected in the nearly equal values of Gin i, Multiplicity, 
andM2o coefficients from the analysis in § |4.2| (see Figure[7j). 
In the details, there are however some noticeable differences. 
In particular, a few of the brighter peaks, or clumps, visible 
in the NIC2 or in the Ha image have no obvious counterparts 
in the other map. The properties of the clumps identified in 
these two images are studied in Paper II. The central regions 
of BX482 also appear to be slightly brighter in H^q band 
emission than in Ha emission. 

For the other galaxies with lower resolution SINFONI data, 
the small-scale substructure is smeared out but the overall 
H\(,o band and Ha distributions at 0."4-0."6 resolution are 
very similar. For BX389, the southern companion that is 
clearly separated from the main galaxy in the original NIC2 
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image is still apparent at lower resolution as an extension in 
the outer isophotes. The well-resolved SE and NW merger 
components of BX 528, and the bridge of emission connect- 
ing them, blend together and the PSF-matched NIC2 image 
strikingly resembles the fairly regular and centrally concen- 
trated Ha line map. While in most cases, the general large- 
scale features of the morphologies are recognizable at pa Of 5 
resolution, BX 528 offers an example where the nature of the 
system can only be reliably assessed with the additional infor- 
mation from either the kinematics or higher spatial resolution 
data. 

5.2. Comparison of Structural Parameters 

We performed 2D Sersic model fits on both the Ha maps 
and the PSF-matched NIC2 ima ges us ing GALFIT, follow- 
ing the procedure described in § 14.1.11 For BX 482, we ran 
fits on the AO-assisted as well as the seeing-limited Ha data. 
For BX528 and BX389, we considered only single compo- 
nent fits since, at the ps Of 5 resolution of the SINFONI data, 
the merger components are unresolved in the former case, and 
the southern companion is blended with the primary galaxy in 
the latter. For consistency, fits on the convolved NIC2 images 
were restricted to a smaller region corresponding to the deep- 
est area with effective Pa 3. "8 x 3. "8 FOV of the SINFONI data 
at the Of 125 pixel" 1 scale, and Pa 2."3 x 2."3 FOV of the AO 
data at Of 05 pixel" 1 for BX 482. This reduced area potentially 
affects the sky background estimates, degenerate with the pro- 
file shape and half-light radius. The smaller FOV as well as 
the broader PSF both contribute to the significantly larger un- 
certainties for the fits performed on the lower resolution maps 
compared to those using the maps at the original NIC2 reso- 
lution. Table [5] reports the best-fit structural parameters and 
Figure QT| compares the results from Ha with those from the 
PSF-matched NIC2 data. 

Overall, there is good agreement between the measure- 
ments based on the different data sets, with differences al- 
most always within the (larger) 1 a uncertainties of the results 
obtained from the lower resolution maps. BX 663 is the most 
deviant source, with largest offsets in Sersic index (smaller by 
1.3— 1.6), half-light radius (smaller by factors of 1.14— 1.35), 
and PA. (differences by pa 15° —25°) between the lower reso- 
lution Ha and PSF-matched H\(& band maps and the original 
resolution H^o band image. On the other hand, the param- 
eters derived from the Ha and PSF-matched //[6o maps are 
in closer agreement, with R e and n for Ha higher by 18% 
and 0.3, respectively, and b/a lower by 13%, all within the 
1 a uncertainties. To investigate the cause of the differences 
in parameters obtained from the original Hj6o band image of 
this galaxy and from the smoothed and resampled version, we 
created a pure n = 2 Sersic model with total magnitude, R e , 
b/a, and PA. equal to the values of the best fit to the original 
NIC2 map of BX663. The mock galaxy was pasted into an 
empty area of the full NIC2 map, which was then convolved 
to the SINFONI PSF, and finally resampled to 0." 125 pixel" 1 . 
The input parameters are recovered within < 1 a by the fitting 
procedure in the Of 05 pixel" 1 data, both at the original NIC2 
and the lower SINFONI resolution (with differences of < 4% 
in R e and < 0.2 in n). Significantly lower R e and n index are 
derived after resampling to the larger pixel scale with, how- 
ever, the best-fit parameters for the mock galaxy in excellent 
agreement with those of the real PSF-matched and resampled 
#160 band image of BX 663 (within < la, with 4% difference 
in R e and 0.12 in n). These tests indicate that the informa- 



tion on the steep inner profile of BX 663 leading to n = 2 is 
lost with the coarser Of 125 pixel" 1 sampling despite the good 
pa 0."4 resolution of its SINFONI data. 

For the other five sources, a comparison of the fits to the 
smoothed and original H\(,o maps indicates that the half-light 
radii, Sersic indices, and P.A.'s agree on average to < 1%, 
0.1, and 4°, respectively, with differences among individual 
galaxies of at most 2%, 0.2, and 7°. The axis ratios are higher 
by pa 30% on average (in the range 5% -55%) for the fits 
to the smoothed NIC2 data. Among these five sources, the 
differences between fits to the pa Of 5 resolution Ha and PSF- 
matched //i6o maps are on average 3% for R e , 0.08 for n, and 
5° for PA. (at most 17%, 0.12, and 11°). The axis ratios from 
the Ha maps are higher than those from the smoothed H\a) 
data by 20% on average (in the range 3%-30%). For BX 482, 
the differences between the parameters obtained based on the 
seeing-limited and AO Ha maps are in the same ranges as 
for the smoothed versus unsmoothed band maps for the 
other sources (excluding BX663). The close agreement be- 
tween the parameters derived from the Ha and Z/i6o band 
emission for BX482, either at pa 0."15 or pa 0f5 resolution, 
is unsurprising in view of their similar morphologies in the 
high resolution SINFONI and NIC2 data (FigureslTlandfTOl). 

One important implication of the Sersic model fits is that the 
global Ha and rest-frame pa 5000 A light distributions do not 
differ significantly from each other. This result is further sup- 
ported by the comparison of the radial light profiles extracted 
along ellipses with center, axis ratio, and PA. from the best- 
fit model to the H\a) band maps in Figure [T2] In this Figure, 
the lower resolution Ha and PSF-matched //i6o maps are used 
for all sources, except for BX 482 where the profiles are ex- 
tracted from the AO-assisted Ha data and original Z/ieo band 
image. At pa 0."5 resolution, the Ha and H\&q band profiles 
follow closely each other, with differences becoming impor- 
tant only beyond sa 1 .5 — 2/? e or typically 8-10 kpc. For the 
low flux levels reached at these large radii, systematic resid- 
uals in background subtraction and flatfielding potentially af- 
fect the profiles in a substantial way. At radii < \.5-2R e , 
some of the (small) discrepancies may reflect differences in 
the detailed morphologies as well as uncertainties in the ex- 
act SINFONI PSF and thus the convolution kernel applied to 
smooth the NIC2 data. For the higher resolution profiles of 
BX482, the ring-like feature that is dominated by the very 
bright clump on the southeast side of the galaxy is promi- 
nent in both Ha and rest-frame pa 5000 A; as seen also from 
Figure [lOl the regions in the vicinity of the galaxy center are 
comparatively fainter in Ha than in //i6o band light. 

The comparisons above do not indicate significant beam- 
smearing effects between angular resolution of FWHM pa 
Of 15 and pa 0f5 for the objects in our NIC2 sample. The 
galaxies are sufficiently large that their overall radial light dis- 
tributions remain well resolved at pa Of 5 resolution. The most 
noticeable difference is the systematic shift to higher axis ra- 
tios for all sources. The coarser Of 125 pixel" 1 sampling leads 
to a significant decrease in Sersic index for BX 663, the object 
with steepest inner light profile (n = 2) at 0."05 pixel" 1 . The 
parametric fits do not reveal any major difference between Ha 
and broad-band rest-frame pa 5000 A emission in terms of R e , 
n, and PA. for the five sources with n < 1. 

5.3. Emission Line Contributions 

In view of the similarity between the Hi^o band and Ha 
maps for our targets, it is important to assess the contribu- 
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tion of emission lines to the broad-band emission. At the red- 
shift of our targets, the NIC2 F160W bandpass includes the 
[O III] AA4959,5007 doublet and H/3. For all sources except 
BX528, SINFONI seeing-limited observations of these lines 
with t he H band grating were take n as part of the SINS sur- 
vey dForster Schreiber et al.ll2009l). These data are analyze d 
and discussed by P. Buschkam p et al. I ( 1201 1. in p reparation), 
to which we refer for details. The [OIII] lines are detected in 
all galaxies but BX663, while H/3 is detected in MD41 and 
BX610. Some of the non-detections can be attributed to the 
galaxies' lines being redshifted to the wavelength of bright 
telluric lines (in particular for BX663 and BX482). For 
all non-detections, we derived upper limits on the fluxes. 

Table [6] reports the fractional emission line contribution 
based on the integrated Hi6o band flux density and line fluxes 
(or 3 a upper limits). In all five sources with SINFONI H- 
band data, H/3 is the weakest line. The largest line contri- 
bution is inferred for BX 389, where the [O III] lines make 
up 24% of the total Hi^) band flux density and, with the 3 a 
upper limit on H/3, the total line contribution is estimated at 
< 29%. BX610 has the smallest contributions, amounting 
to rj 6% for all three lines combined. The broad-band H\a) 
emission of our NIC2 targets is thus clearly dominated by the 
continuum, with only a modest to very small fraction from the 
[O III] and H/3 lines. 

While the general appearance of the targets in [O III] emis- 
sion, and in H/3 where detected, resembles that in Ha, there 
are in some cases measurable differences. These gradi- 
ents in line ratios reflect spatial variations in line excita- 
tion or gas-phase o xygen abundance, and in dust extinction 
across the targets dGenzel et al. 1 12008; P. Buschkam p et al.l 
1201 1. in preparationl) . To assess possible spatial variations in 
the emission line contributions, we also examined the galax- 
ies on a pixel by pixel basis after convolving the NIC2 maps 
to the PSFs of the SINFONI H band data, resampling them 
to 0."125 pixel" 1 , and aligning them following the same pro- 
cedure as for Ha. We considered only pixels with S/N > 3 
in both //i6o flux density and line flux (the summed fluxes 
over these pixels represent typically > 75% of the integrated 
measurements of the galaxies in 3 "-diameter apertures). The 
mean and median of the [O III] A 5007 contributions for the 
four sources where the line is detected are within w 7% of the 
estimates from the integrated fluxes. The top (75%) quartile of 
the pixel line contributions range from 7% (BX610) to 16% 
(BX389). For H/3, detected in MD41 and BX610, the mean 
and median contributions in individual pixels are within < 5% 
of those determined from the total fluxes, with top quartiles of 
7% for MD41 and 2% for BX610. We conclude that spatial 
variations in line emission are not likely to affect significantly 
the observed H\a) band morphologies. 

6. COMPARISON WITH OTHER Z ~ 2 SAMPLES 

As we found in §|4] there are hardly any significant trends of 
structural parameters or morphological properties with global 
stellar and dust extinction properties among our six SINS 
galaxies. However, the small sample size leads to limited 
statistics and range in parameter space covered. In order 
to expand our analysis and ultimately place the morpholog- 
ical results for our NIC2 sample of rest-UV selected galax- 
ies in the larger context of the z ~ 2 galaxy population, we 
compare their properties to those of galaxy samples selected 
using very different photometric criteria. We restrict this 
comparison to samples in the same redshift range that have 
been characterized using the same type of observations, i.e. 



HST/NIC2 F160W imaging. The two c omparison samples 
are 19 A"-s elected galaxies pr esented by Ivan Dokkum et al.l 
(2008) and iKriek etaLl ( 120091). and 11 24 /xm-s elected dusty, 
IR-luminous galaxies from lDasvra et al. I (120081). 

The paren t samp le for the Ivan Dokkum et alH (2008) and 
IKriek et aD (120091) analyses is described by Kriek et al. 
(120061 12008a) and consists of 36 massive galaxies with 
^vega < 19.7 mag, 1.6 < z < 3.0 and log(MjM^) > 10.5. 
The stel l ar ma sses are given by Kriek et al. (2008a) for a 
Salpeterl (119551) I MF; we d i vided them by 1.7 for consistency 
with our adopted Chabrier (2003) IMF throughout the analy- 
sis below. All galaxies in this /f-selected sample have spec- 
troscopy with the Gemini near-IR spectrograph (GNIRS), and 
19 objects obtained early in the survey, at 2 < z < 2.6, have 
NIC2 imaging with the F160W filter. Based on their GNIRS 
spectra, nine of these objects present no detectable Ha emis- 
sion, with an upper limit of VK(Ha) = 10 A. The remaining 
ten objects have detectable rest-frame optical emission lines, 
indicating either star form ation or AGN activ ity (for six and 
four objects, respectively: IKriek et al.l 120091) . The range of 
activity in this roughly stellar mass-selected sample of 19 ob- 
jects is reflected in their structural parameters. The nine qui- 
escent galaxies are very compact, nearly circularly symmet- 
ric, and show little if any substructure except for one source. 
The y a re characterized by a median circularized effective ra- 
dius™ of 0.9 kpc (ranging from 0.5 to 2.4 kpc), and a me- 
dian Sersic index, n, of 2.3 (ranging from 0.5 to 4.5). The 
six star-forming galaxies are signicantly more extended and 
less centrally-concentrated, and appear more similar to our 
SINS targets. Their median R e . c hc is 2.8 kpc (ranging from 
1.5 to 5.0 kpc) and median n index is 0.9 (ranging from 0.0 
to 3.3). The four AGN sources resemble the quiescent galax- 
ies on average in terms of structure, with a median size of 
^e.circ ~ 1-1 kpc and Sersic index of n ~ 2.5. 

As discussed by Krie k et al. I d2009), in terms of estab- 
lishing evolutionary connections within this /f-selected NIC2 
sample, half of the star-forming galaxies appear to consti- 
tute direct structural precursors to the AGN and quiescent 
galaxies, based on the fact that the majority of their stellar 
mass may be contained within a compact core. However, the 
remaining star-forming galaxies lack evidence for a strong, 
central mass concentration and t herefore any direc t structural 
link with the quiescent objects (IKriek et al.l 120091) . The me- 
dian stellar mass for the quiescent galaxies is 1.6 x 10 11 M Q , 
implying central stellar mass densities ^100 times higher 
than for z ~ galaxies with similar stellar mass, character- 
ized by significantly larger effective radii (~ 5 kpc). Under- 
standing the formation and evolution of such compact, qui- 
escent galaxies at high z is the basis o f much ongoing ob- 
servational and the oretical work (e.g.. iTruiillo et alH 120071; 
Cimat ti et al.ll2008l;lyan der Wei et al ll2008HDamianov et al l 
2009 1 ; iHopkins et al.ll2009tlMancini et al.ll2010tlWuyts et al.l 
2010). 

The paren t samp le for the objects we analyzed from 
iDasvra eta l. (2008) consists of 52 IR-bright sources in the 
Spitzer Extragalactic First Look Survey, with 24 /im fluxes 
greater than 0.9 mJy, i//v(24^m)/^/v(8/im) > 0.5, and 
^/v(24/im)/^/v(6440A) > 1. The latter two criteria were 
designed to select dusty starbursts at z > 1, or obscured 
AGNs. Thirty-three such sources at z > 1.5 were imaged 
with NIC2 and the F160W filter, 11 of which are in the 

14 The circularized effective radius is defined as ff ec ire = yb/a. 
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2.0 < z < 2.5 range of interest here. All of these sources 
have rest-frame 14 /im luminosities greater than 10 12 L Q (de- 
rived from Spitzer/IRS data), implying total IR luminosities 
greater than ~ 10 13 L Q . Each is classified as one of mul- 
tiple types based on its mid-IR IRS spectrum. The sample 
considered here contains the following types: "PAH," indi- 
cating a spectrum dominated by PAH emission (one object); 
"AGN," indicating a spectrum characterized primarily by the 
AGN continuum (two objects); "obscured," indicating a spec- 
trum dominated by the SiO absorption feature at 9.7 /im (four 
objects); or, finally, "mixed," in which no single starburst or 
AGN component dominates the spectrum (four objects). The 
median NIC2 H\§o band magnitude for these 11 objects is 
#160,ab = 22.19 mag, ranging from 21.51 to 24.10 mag, and 
their median half- light radius (estimated u sing the SExtractor 
software package; iBertin & Ar nouts 1996) is 2.4 kpc, ranging 
from 1 .5 to 4.3 kpc. These sources show varied morphologies, 
from compact and regular to multiple systems, to disturbed 
and/or irregular; the majority, however, tends to be more com- 
pact a nd concentrated than our SINS targets. iDasyra et al.l 
(2008) assign a qualitative description of the morphology of 
each object in their NIC2 sample, to assess the importance 
of merger events. Of the 1 1 objects described here, five are 
classified as exhibiting no signs of interaction, five are "Dis- 
torted," and one is "Binary." 

These comparison samples observed with the same NIC- 
MOS camera and filter, and spanning the same redshift range 
as our SINS NIC2 targets, enable in principle a controlled 
comparison of morphological properties. However, slightly 
different analysis techniques, integration times, and pixel 
sampling are featured in the published work. These differ- 
ences have consequences for the derivation of the morpho- 
logical properties, which depend on the depth of the data as 
well as on the pixel sampling and S/N threshold applied for 
the analysis. We therefore obtained the actual images for 
these samples (van Dokkum 2008, private communication; 
Yan 2008, private communication; Kriek 2010, private com- 
munication) and performed a uniform analysis of our NIC2 
sample alongside each of the comparison samples. 

6.1. Comparison of Structural Parameters 

As a first step in our comparison, we considered the global 
structural parameters and, more specifically, the effective ra- 
dius R e and Sersic index n. While each of the SINS targets 
was observed for four orbits, 17 of the 19 ^-selected tar- 
gets were observed for three orbits and the two brightest ob- 
jects for two orbits. In the 24 /im-selected sample, five ob- 
jects were observed for two orbits, and s ix for only one or- 
bit. M oreo ver, the publ i shed i mages from lvan Dokkum et al.l 
(2008) and lKriek et all (120091) wer e drizzled to 0."038 pixels, 
while those of lDasyra et al. I (120081) to 0."076 pixels. To match 
both the imaging depth and sampling of the two comparison 
samples, we created 1-, 2-, and 3-orbit subsets of the NIC2 
data of our SINS targets, and also drizzled them to 0."038 
and 0."076 pixel" 1 . We repeate d our s ingle-component Sersic 
model fitting as described in § 14. 1 . 1 I for each combination of 
exposure time and pixel sampling. Because of the differences 
in d epth and pix el sampling of the Ivan Dokkum et al.l (2008) 
and IKriek et al.l (120091) NIC2 data with respect to those of 
IDasyra et al. I d2008l) . we compare our SINS targets separately 
with the ^-selected and 24 /x m-selected samples. 

For the .g - select ed galaxies. Ivan D okkum et al. (200§1) and 
IKriek et aTl (120091) derived the structural parameters using the 



same code as we employed (GALFIT) and a conceptually 
analogous approach was followed in treating the sky back- 
ground. We thus adopted the best-fit R e and n parameters 
as published by these authors, converting the circularized ef- 
fective radii to standard effective radii. None of our conclu- 
sions is affected if we convert our radii to circularized radii 
instead. For the SINS galaxies, the difference in structural 
parameters between the 3- and 4-orbit subsets is negligible. 
The pixel scale has fairly little impact on the derived pa- 
rameters of extended, well-resolved sources but more so for 
compact ones. Accordingly, for all but one of our SINS tar- 
gets, the resulting best-fit R e and n obtained from the 0."038 
and 0."05 pixel" 1 images are nearly identical and within the 
(small) 1 - 2 a uncertainties, depending on the source and pa- 
rameter. For BX663, the AGN host with a compact central 
peak, the changes are significantly larger than for the other 
five SINS sources although they are formally within 1 .5 a of 
the best-fit values for each pixel scalePl 

For the 2 4 /xm-s elected dusty IR-luminous sample, 
IDasyra et al. I (120081) did not perform single-component 
Sersic profile fitting for their analysis. We thus applied our 
fitting procedure using their images to obtain the R e and n 
for comparison with our sample. The resulting effective radii 
agree on average within 20%, or < la based on our fitting 
uncertainties , with those estimated using the SExtractor 
software by IDasyra et al. I . The most important difference 
(by 30%, or 3 a) occurs for the largest, most diffuse and 
asymmetric source (MIPS 8242). For our SINS targets, and 
compared to the nominal fits for the 4-orbit 0."05 pixel" 1 data, 
the differences resulting from the shallower subsets and larger 
0."076 pixel" 1 sampling are typically < 3 a. Again, BX663 
shows the largest changes with R e smaller by « 0.6- 1.2 kpc 
and n lower by sa 0.7 - 1, depending on the subset. 

Figure Q~3] plots the best-fit n index versus R e of our SINS 
sample along with those of the /^-selected and 24 /im-selected 
comparison samples. In terms of size and Ser sic index, 
the quiescent sample of Ivan Dokkum et al.l (120081) is clearly 
distinct from our SINS galaxies, with overall significantly 
smalle r R e and typically higher n values. The Krie k et al.l 
(2009) objects are in termediate between the passive objects 
of Ivan Dokkum et al. I and our acti vely star-for ming SINS 
sources with, as already discussed bv lKriek et al. L a clear sep- 
aration between the AGN sources that have sizes and n in- 
dices in a narrow range and fully overlapping with the passive 
objects whereas the purely star-forming sources span wider 
ranges shifted towards larger sizes and lower n values. 

Quantitatively, our SINS objects lie at the large size - low 
n end with sample median /? e ,med = 5.1 kpc and « me d = 0.35 
(or 5.7 kpc and 0.33, respectively, when excluding the AGN 
source BX 663). In the ^-selected comparison sample, the six 
star-forming sources have median values of /? e .med = 4.3 kpc 
and n me d = 0.95, the four AGN sources have ^ e .med =1-4 kpc 
and « me d = 2.46, and the nine passive objects have ^ ejme d = 
1.5 kpc and n me d = 2.30. Our AGN source BX663 is not 
as compact as the IKriek et aTl AGNs but, among our SINS 
sample, it reflects the trend towards lower R e and higher n of 
AGN relative to star-forming galaxies see n in the .g-selected 
sample. One source (HDFS1-1849) from lvan Dok kum et al. 
appears to be an "outlier" with respect to the other eight pas- 

15 Specifically, the best-fit i? c and n for the finer 0."038 pixel -1 sam- 
pling are 3.6;^ 3 kpc and 1.31^?'J|, respectively, compared to | kpc and 
^ for the 0"05 pixel -1 one. 
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sive ^-selected objects: its R e = 4.4 kpc and n = 0.5 is c learly 
more characteristic of the SINS and iRriek et al. I ( 120091) star- 
forming galaxies. As we will see below, it is also more simi- 
lar to the star-forming sources in terms of the non-parametric 
Gini, and M20 coefficients. This object is anomalous 
within the sample of quiescent galaxies insofar as its broad- 
band SED fit indicates in fact active star formation. The lack 
of emission lines that formed the basis of the "quiescent" clas- 
sificati on may only be d ue to this object's dusty, edge-on na- 
ture (IRriek et al. 1 120091) . Excluding this object, the median 
values for the eight remaining passive /f-selected sources are 
^e.med = 1 . 4 kpc and w = 2.55. 

With the Dasv ra et al.l (2008) sample, we can explore vari- 
ations in different regimes of bolometric luminosity. While 
the four SINS objects with Spitzer/MlPS data are among the 
most luminous of the rest-UV selected high z population, their 
observed 24 pm fluxes and total IR luminosities are still on 
average almost an order of magnitude lower than those of 
the 24 /im-selected sources of iDasvra et all (120081) . Overall, 
the 24 //m-selected dusty IR-luminous sample is distinct from 
our SINS objects, and appears to be intermediate between the 
rest-UV SINS and ^-selected star-forming galaxies and the 
^-selected passive and AGN objects. The 24 /im-selected 
sources have a median size and Sersic index of R s>me a = 
2.6 kpc and w me d = 1.45. These values are comparable to the 
median rest-frame optical structural parameters derived from 
NIC2 imagi ng for a sample o f high -z submillimeter-selected 
galaxies by Swinba nk et al. | (|2010h. wh ich overlap in bolo- 
metric luminosity with the IDasvra et al. I objects. Sampled at 
the same 0"076 pixel" 1 scale, and 1- and 2-orbit integration 
times, the SINS sample has about twice larger /? e ,med = 5.4 kpc 
and significantly lower n me d = 0.37. The trend for the 24 /im- 
selected sources to lie towards smaller radii and steeper in- 
ner profiles with respect to the SINS ga laxies may refl ect the 
prevalence of AGN activity among the Das yra et al. I sample 
(10 out of 11 objects). 

6.2. Comparison of Morphological Coefficients 

Clearly, the various samples probe different regimes in 
structural parameters, and this appears to be related to both 
the level of star formation and AGN activity. To further quan- 
tify the differences between our SINS galaxies and the com- 
parison samples, we analyzed the image s through the non- 
parametric approach described in § 14.2.11 using the relevant 
data sets for the SINS targets matching the depth and pixel 
sampling of the K- and 24 /im-selected samples. In addition, 
to take into account the redshift range our targets and remove 
the effects of surface brightness dimming, particularly impor- 
tant in this approach, we selected pixels for our analysis as if 
all objects were observed at the maximum redshift among the 
three samples, i.e. z = 2.6. 

We estimated Gini, '5, and M20 for galaxies in the three 
samples, and the results are shown in Figure [14] All panels 
contain our SINS galaxies, and the top two panels show G 
versus \& and G versus M20, using the /f-selected sample for 
comparison. The bottom two panels show the same param- 
eters, but feature the 24 /im-selected sample for comparison. 
For the top two panels, the SINS targets have been reanalyzed 
using 3- and 2- orbit subsets of the data at Of '03 8 pixel" 1 , in or- 
der to matc h the depth and samp ling of the lvan Dokkum et al.l 
(2008) and iKriek et al.l {2009) images. Given the compact 
and axisymmetr ic nature of the qui escent /f-selected subsam- 
ple featured in Ivan Dokkum et aU, we applied the same non- 



parametric analysis to our em piric al PSF, constructed from 
multiple stars as described in § !2.5l but from our 0."038 pixel 
set of images. This analysis indicates that while the quies- 
cent, ^-selected objects are extremely compact, the majority 
of them can be distinguished from unresolved point sources, 
with significantly lower Gini values. For the bottom two pan- 
els, our SINS targets have been reanalyzed using 2- and 1- 
orbit subsets of the data at 0." 076 pixel" 1 , in order to match 
the depth and sampling of the IDasvra et al. I (120081) images. 
The non-parametric statistics were calculated as well for the 
empirical PSF sampled with this pixel scale. In all panels, 
the shallower 3-, 2-, and 1 -orbit SINS data subsets result in 
lower Gini coefficient for a given object, while higher values 
are obtained for '5 and M2o- 

The non-parametric morphological analysis reinforces the 
qualitative visual impression from the images and the re- 
sults from the Sersic profile fitting of the different samples 
of galaxies. The Zf-selected passive and AGN sources oc- 
cupy a disjoint region of G— \f r — M20 morphological parameter 
space from that of the SINS targets. They are virtually all off- 
set towards higher Gini, and both lower if? and M20 values, 
with sample median values of G me d = 0.54, ^? me d = 0.8, and 
A^20,med = _ l-8 for the passive galaxies and essentially identi- 
cal G med = 0.53, * me d = 0.8, and M 20 ,med = -1-8 for the AGN 
sources. As was the case for the structural parameters, only 
a single galaxy from the Zf-selected passive sample (HDFS1- 
1849) even approaches the lower Gini (median of 0.28), and 
higher '5 (median of 7.8) and M20 (median of -0.5) values 
of our SINS sampleQ In contrast to the passive and AGN 
objects, the six /f-selected star-forming objects show consid- 
erably more overlap with the SINS sample. The median val- 
ues for the /^-selected star-forming sources are G me d = 0.34, 
^med = 3.6, and M2o,med = _ l-4, which are intermediate be- 
tween the parameters describing the SINS sample and the K- 
selected passive and AGN objects. 

One slight caveat to these results is the difference in typ- 
ical S/N among the samples, with the ^-selected quiescent 
and AGN objects ~ 2 - 3 times higher in average surface 
brightness than the SINS targets, and the Zf-selected star- 
forming sources ~ 1. 5 times brighte r . As discussed by 
iPeter et al.l (120071) and lAbraham et al.l (120071) . using a fixed 
surface brightness threshold (effectively what we have done), 
will lead to a "truncation bias" in the Gini coefficients esti- 
mated for fainter systems, in that their pixel flux distribution 
functions will be sampled over a smaller dynamic range than 
those of more luminous objects. To test for the significance of 
this effect, we applied a higher surface brightness threshold to 
the ^-selected samples, such that the resulting pixel flux dis- 
tribution functions span the same dynamic range as those of 
our SINS targets. We find G me d = 0.45 for the truncated pixel 
flux distribution functions for Zf-selected passive and AGN 
sources, still significantly higher than the median of the SINS 
sample, and still with no overlap in morphological parameter 
space. Applying the same type of truncation to the A'-selected 
star-forming galaxies, we find very little effect on the sample 
median properties, with G me d = 0.33. 

Turning to the properties as a function of bolometric lu- 
minosity, in the G — \& — M20 parameter space, the 24 /im- 
selected dusty IR-luminous sample separates on average from 
the SINS sample in the same sense as the /^-selected sample, 
if not as significantly as the passive and AGN objects. The 

16 These median values repre sent the ensemble of 3- and 2-orbit subsets 
of the SINS data shown in Figure [l4l 
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median parameter values for the 24 /im-selected dusty IR- 
luminous sources are G me d = 0.44, \IVed = 2.1, and M2o,med = 
—1.2. Sampled at the same pixel scale and exposure time, the 
SINS targets have lower Gini (median of 0.30), and higher 
(median of 6.4) and M20 (median of -0.6). The typical S/N 
of objects in the two samples is very similar, so no issue of 
truncation bias applies in this comparison. These quantita- 
tive morphological results confirm again the qualitative visual 
impression that our SINS targets have systematically more ir- 
regular, extended, and clumpy morphologies tha n those of the 
24 /im -selected dusty IR-luminous sample of Dasv ra et al.l 
(200§. 

Luminous dusty systems with Lir > 10 12 L are com- 
monly associated with m ajor mergers in the local universe 
(e.g.. iGenzelet al.l l2Q0Th . The 24 /im-selected sample of 
Dasyra et al. (2008) falls well above this luminosity thresh- 
old, yet the processes giving rise to extreme bolometric 
luminosities at high redshift may be more div erse, based 
on theoretical grounds (e.g., iDave et all 120101), and sup- 
ported by observation al evidence (e.g ., [Dasyra etaDl2008l: 
ISwinbank efaTI 120101) . iDasvra et al.l find roughly 50% of 
their full NIC2 sample to appear morphologic ally "quies- 
cent." Furthermore, consistent with the results of Rrie k et al.l 
(2009), these authors show that objects in their sample with 
mid-IR AGN signatures also tend to have more compact 
rest-frame optical morphologies than those with starburst- 
dominated spectra. Given that all but one object in the 
IDasvra et al.l z = 2.0-2.5 sub-sample show evidence of AGN 
activity, the 24 /im-selected dusty IR-luminous galaxies in- 
cluded in Figure [T4l are among t he more com pact objects 
in the full Dasvra et al. sample. Dasv ra et al. I hypothesize 
that 24 /im-selected IR-luminous objects, with warmer mid- 
IR/far-IR colors th an the submillimeter-s elected IR-luminous 
high z population (Chap man et al.|[2003l) . may represent the 
merged, dynamically-relaxed counterparts of the (major-) 
merging submillimeter galaxies. On the other hand, the 
rest-frame optical structural param eters and sizes of submil- 
limeter galaxies as presented by Swin bank et al. I d2010b are 
very similar to those of the 24 /im-selected IR-luminous ob- 
jects. A more careful morphological comparison of the two 
bolometrically-luminous samples is clearly required to con- 
struct an evolutionary sequence. 

6.3. The Connection Between Structure and Stellar 
Populations in z ~ 2 Galaxies 

In analogy with trends seen in the local universe 
dKauffmannet al.l 120031: lShen etal.1 [2003), the structural 
properties of z ~ 2 galaxies appear to be correlated with 
their stellar populations. B ased on imaging with both 
NICM OS/NIC3 from space (IZirm et al. I 120071; i Toft et al. 
20071) . and with VLT/ISAAC dFranx et al.ll2008t Toftetal. 
2009) and UKIRT/WFCAM (Williams et al. 1120 lob from the 
ground, it has been shown that quiescent, high redshift galax- 
ies have systematically smaller sizes and higher stellar mass 
surface densities than actively star-forming systems. Based 
on the sample of 1 9 ZT-selected galaxi es discussed in the pre- 
vious subsections, iRriek et ail d2009l) used their higher res- 
olution NICMOS/NIC2 imaging to demonstrate that specific 
SFR and rest-frame U — B color are strongly correlated with 
galaxy structural parameters. We here examine the struc- 
tural properties in connection with the ste llar populations for 
our SI NS s ample, complemen ted with the Ivan Dokkum et alH 
(2008) and iKriek etaU (120091) ^-selected samples for which 
stellar properties from SED modeling are available. 



As proxies for galaxy stellar populations, we consider the 
stellar mass and specific SFR. Within our SINS NIC2 sample, 
the absence of correlation between morphology and stellar 
population properties includes parameters derived from Sersic 
models (R e and n) as well as the non-parametric morphologi- 
cal coeffici ents (G, jfr , and M 2Q). This result echoes previous 
ones from lLawet al.l d2007al) and lPeter etaTI (12007b demon- 
strating a lack of correlation between stellar populations and 
rest-UV morphologies f or UV-selec t ed ga laxies at z ~ 2. 
However, as described by Pet er et al. I d2007) the UV-selected 
galaxies do not span the full range of stell ar population prop- 
erties that exist at z ~ 2 (see also, e.g.. Ivan Dokkum et al.l 
l2006tlM. Damen et al. 11201 1. in preparation!). The g-selected 
sam ple with NIC2 im aging from Ivan Dokkum et all (120081) 
and lKriek etaD (120091) ideally augments our SINS NIC2 sam- 
ple in this respect, considerably extending the range probed 
in specific SFR by over two orders of magnitude, in the same 
redshift range and at comparable stellar masses as our SINS 
targets. 

Figures [15] and [16] combine the morphological properties 
of our SINS and /f-selected samples, now plotted as a func- 
tion of the stellar mass and specific SFR derived from SED 
modeling. The SINS and /f-selected samples overlap in stel- 
lar mass over the range from w 3 x 10 1() to s» 10 11 M Q . Over 
this range, there is significant diversity in structural parame- 
ters, with the SINS targets exhibiting the largest R e , smallest 
n, largest and M20, and smallest G values. At the other ex- 
treme, we find the AGN and quiescent /f-selected galaxies, 
with the star-forming ^-selected galaxies showing intermedi- 
ate properties. More specifically, over this interval in stellar 
mass, the galaxies cover a wide range in R e , consistent with 
the spread of about an order of magnitude in rest-frame opti- 
cal sizes for z ~ 2 galaxies of sim i lar masses found in previous 
studies (e.g . iTruiilloet al.ll2006t IZirm et al.ll2007l: lToftetal.1 
20071 120091: iFranxet al.ll2008l) . 

On the other hand, the /^-selected and SINS samples sep- 
arate more cleanly in specific SFR. A clear sequence of in- 
creasing specific SFR is represented, respectively, by the qui- 
escent and AGN ^-selected galaxies, the star-forming K- 
selected galaxies, and the SINS sample. In turn, signifi- 
cant connections with specific SFR are found for every mor- 
phological coefficient considered here, with systematically 
higher R e , "f, and M20, and systematically lower n and G, 
as the specific SFR increases. These strong trends with 
star formation activity are c o nsistent with the p r evious work 
fromlToft et al.l d2007Ll2009l) ; IZirm et all d2007l) : lFranx et al.l 

(2008) , but are based on higher resolution NIC2 imaging. 
Furthermore, they extend towards higher specific SFR values 
the NIC2-based c orrelations between morph olog y and stellar 
popula tions from Ivan Dokkum et alH d2008l) and lKriek et al.l 

(2009) . The re s ults p resented here reinforce the conclusion 
by iFranx et al.l d2008l) that galaxy structural parameters are 
better correlated with specific SFR and color than with stellar 
mass. These results highlight the fact that stellar mass alone is 
not a good predictor of structural parameters, nor of detailed 
morphologies as we find from our non-parametric morpho- 
logical analysis. 

In that respect, it is also interesting to consider morpholo- 
gies as a function of the stellar mass surface density of the 
galaxies. Indeed, in the local universe, the specific SFR (as 
well as rest-frame optical colors) is observed to be much more 
closely rel ated to stellar mass surface de nsity than to stellar 
mass (e.g.. lKauffmann et al.ll200l 120061 for the SDSS sam- 
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pie), a nd these relationships are found to pers ist out to at least 
Z ~ 2 dFranxet all 12008: Maier_etaL][2009). This trend is 
again reflected among the SINS and ^-selected samples dis- 
cussed here. Figure [17] plots the specific SFR as a function 
of the projected stellar mass surface density of the ob- 

jects. We calculated the surface density assuming that the ob- 
served 1 .6 light (rest-frame 5000 A) from the NIC2 im- 
ages traces well the stellar mass (as we will demonstrate for 
one of our SINS targets in § 17. It , and thus that half of the to- 
tal stellar mass is enclosed at r < R e . We find a clear trend of 
decreasing specific SFR with increasing mass surface density, 
which is fu lly consisten t with the median relationship at z ~ 2 
derived by Fra nx et aTl d2008). Our six actively star-forming 
SINS targets lie at the lower end of projected mass surface 
density with median of log(£(Af*) [M kpc~ 2 ]) = 9.0 (or 8.8 
when excluding the AGN source BX 663) a nd the nine quies- 
cent galaxies of Ivan Dokkum et al. I (120081) lie at the higher 
end with median of log(S(M*)[M Q kpc- 2 ]) = 10.5 (or 10.6 
when omitting the a nomalous HDFS1- 1849). The ten emis- 
sion line galaxies of Krie k et al. I (|2009) span the intermediate 
range, with sample median of log(S(M i )[M Q kpc" 2 ]) = 9.5. 
Among these objects, the subsamples of six star-forming and 
four AGN sources have, respectively, lower and higher mass 
surface densities (median values of log(S(M*) [M©kpc~ 2 ]) = 
9.2 and 10.3). Since the structural and morphological proper- 
ties of the SINS and /f-selected samples considered here are 
strongly connected with the specific SFR, it follows from the 
specific SFR versus £(M*) trend that they are also closely re- 
lated to the stellar mass surface density. We note here that 
the hig h-z submillimeter-selected galaxies of lSwin bank et al. 
(2010) have properties implying on average both higher spe- 
cific SFR and £(M*) than the S INS galaxie s , and d o not fol- 
low the median relationship of iFranx et al. I (12008b . pointing 
towards differences in their assembly histories. 

7. DISCUSSION 

7.1. Rest- frame UV Structure, Colors, and Mass Distribution 

ofMD41 

In this section, we focus on MD41, the target among our 
NIC2 sample with additional high resolution optical imaging 
with HSTI ACS. We first combine the ACS and NIC2 images 
to compare the structural parameters derived from the rest- 
frame UV and optical surface brightness distributions. We 
then use the observed optical-to-near-IR color map to con- 
struct the spatially-resolved stellar mass map of the galaxy. 

7.1.1. Rest-frame UV Structural Parameters 

In § 14.21 we presented the ACS F814W image that exists 
for one of our NIC2 targets, MD41. The fairly similar ap- 
pearance of this galaxy in the ACS igu and NIC2 Hi6o bands 
(Figure [HJ is supported by our quantitative non-parametric 
analysis, which yields very comparable G, VP, and M20 val- 
ues for both images. Here we extend the comparison to the 
structural parameters between the H160 and ig u ban d images. 
This comparison also complements the one in § 15.21 where we 
found that the size, Sersic index, and axis ratio are very simi- 
lar between the Ha and rest-frame m 5000 A continuum emis- 
sion for all our targets at a common (PSF-matched) w 0"5 
resolution. At the redshift of MD41 (z = 2.17), the ACS 
i8i4 band probes rest-frame w 2600 A light, predominantly 
tracing young massive stars similarly to Ha. Therefore, the 
ACS+NIC2 imaging allows us to examine the relative surface 



brightness distributions of different stellar populations at sig- 
nificantly higher resolution than possible with the SINFONI 
seeing-limited Ha data available for MD41, although with 
the caveat that extinction effects are more important at rest- 
frame sa 2600 A than at the ~ 5000 -6500 A mapped by the 

#160 band and Ha observations. 

We repeated our Sersic model fitting procedure (§ 14.1. It on 
the ACS 2gi4 band image of MD41. For consistent compari- 
son with the results from the NIC2 H160 band image, we first 
convolved the ACS data to the spatial resolution of the NIC2 
data. We constructed an empirical PSF for the ACS data using 
ten bright but unsaturated, isolated stars identified across the 
full map. From a Gaussian fit, the PSF FWHM is 0."099 with 
ellipticity of 0.02 (or FWHM = 0"080 and same ellipticity 
for a Moffat profile fit). The smoothing kernel was computed 
through a Lucy-Richardson deconvolution algorithm, and the 
curves-of-growth of the NIC2 and smoothed ACS PSFs agree 
to better than 10% at radii 0."1 < r < l'.'5. No spatial resam- 
pling of the data is needed since both ACS and NIC2 data 
sets were drizzled to 0."05 pixels. For accurate registration 
between the NIC2 and ACS maps, we relied on the star and 
a number of very compact sources visible in both and 
i8i4 band emission within the smaller NIC2 FOV. The right 
panel of Figure [8] shows an RGB color composite with the 
//160 band image in red, and the PSF-matched /814 band image 
in green+blue. 

The best-fit parameters based on the smoothed ACS map 
are: R e (hu) = 7.44jft» kpc, n(ig U ) = 0.27!°;*j, b/a(ig U ) = 
0.30^3:02' and P-A.(/gi4) = 38.8!} i deg. Compared to the pa- 
rameters derived from the NIC2 image (Table|4|, the effective 
radius is thus ss 30% larger, the Sersic index is m 0.25 lower, 
and the axis ratio and PA. are virtually identical. The larger 
size and lower light concentration index derived from the i'sn- 
band image could be due to higher extinction in the central 
regions and/or more extended distribution of the young stel- 
lar population dominating the rest-UV emission. While the 
differences in derived size and n index are significant with re- 
spect to the uncertainties, the fits to the 2314 and band 
images are in broad agreement. Both indicate that MD41 
is a large and fairly edge-on system (implied disk inclina- 
tion angle of m 72°), with a shallow inner light profile. The 
structural analysis supports the conclusions from the non- 
parametric morphological analysis that there is no strong k- 
correction between the rest-frame UV and optical morpholo- 
gies for MD 4 1 . This result is consistent with the general find- 
ings for rest-UV selected samples at z ~ 2, and in contrast to 
those for samples selected at longer wavelengths with typi- 
cally redder rest-frame optical colors ([Dickinson et al."ll2003t 
[ Peter et al.ll2007HLaw et al. Il2007at lToft et al. 1120071 see also 
lOverzier et al. 1 12010 for rest-UV selected samples at z ~ 3-4 
and rest-optical selected samples at z ~ 2). 

7.1.2. Rest-frame Colors, Mass-to-Light Ratio, and Mass 
Distribution 

We can further exploit the ACS and NIC2 images to inves- 
tigate spatial variations in observed colors and mass-to-light 
ratio, and derive the stellar mass distribution within MD 41 . In 
Figure[l8] the left panel shows the observed igu— Hi^o color 
map, over the region for which the pixel flux densities have 
S/N > 3 in each of the ACS and NIC2 map. The observed 
colors are fairly uniform across most of the disk of MD41, 
with mean and median of 0'814 - #i6o)ab = 0.8 and 0.9 mag, 
respectively, and standard deviation of 0.37 mag. This com- 
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pares reasonably well with the integrated color of 0.63 mea- 
sured in a circular aperture of 3" diameterFT. Only two small 
regions exhibit significantly different colors: the southwestern 
edge of MD41 at the location of the two brightest clumps in 
z'814 emission with bluer colors (0'814-#16o)ab ~ 0.2 on aver- 
age), and an area around the center of the galaxy characterized 
by redder colors ((^14— -Hi6o)ab ~ 1.2 mag on average). 

To interpret the color map of MD41, we turn to synthetic 
colors and M/L ratios calcula ted for the z = 2.2 of MD41 
from iBruzual & Charlotl (120031) models for a range of ages, 
extinction, star formation histories (SFHs), and two stellar 
metallicities (solar and 1/5 solar). The details are given in 
Appendix lAl Stellar age, extinction, and SFH are highly de- 
generate in observed iu4~H\60 colors (or rest-frame 2600 A 
- 5000 A colors at the redshift of MD41). However, there 
is a fairly well defined relationship between the stellar mass 
to rest-frame optical light ratio (with luminosity uncorrected 
for extinction) and the observed colors, with very little depen- 
dence on our choice of metallicity. For our purposes, we used 
model predictions through the SDSS g-band filter, a good ap- 
proximation to the rest-frame pa 5000 A luminosity probed by 
the NIC2 F160W filter (see AppendixE). 

We thus derived a mean relationship between the ob- 
served (uncorrected for extinction) M + /L r ? est ratio and ob- 
server's frame igu~H\eo colors for z = 2.2. We used this 
relationship to compute the M*/L* est ratio map of MD41, 
shown in the middle panel of Figure [18] Unsurprisingly, 
the M+/L* est ratio map looks similar to the 1*814— #160 map, 
given the roughly linear relation (in logarithmic units for the 
M/L). The mean and median values of individual pixels are 
log(M*/L" st [MoLg g]) pa -0.9, with standard deviation of 
pa 0.2 dex, The observed integrated color gives nearly the 
same value, -0.83, while the model with best-fit parameters 
from the optical/near-IR SED modeling for MD41 (Tabled 
has a somewhat lower value of -1.08. An average value of 
ps— 1.3 is measured in the regions of bluest colors at the south- 
western edge of MD41, and of pa -0.7 in the redder regions 
around the center. The fairly small scatter among individual 
pixels implies roughly constant observed M*/L r ? est ratio across 
most of the galaxy. 

Multiplying the observed M+/L* est map by the //)6o band 
image as a proxy for the extrinsic rest-frame g-band luminos- 
ity then leads to the spatial distribution of stellar mass across 
MD41. The map of projected stellar mass surface density 
£(M*) is shown in the right panel of Figure [18] It is re- 
markably similar to the //[6o band image itself, indicating that 
the observed 1 .6 /im emission does trace closely the stellar 
mass distribution in MD41. The typical (mean, median) val- 
ues of individual pixels is log(E(M + [M Q kpc" 2 ])) pa 8.0, with 
standard deviation of 0.3 dex. At the largest radii probed, 
and especially on the southwestern edge, the surface den- 
sity drops to log(S(M i )) sa 7.5, and it increases towards the 
center by about an order of magnitude to log(£(M*)) pa 8.4. 
The deprojected surface density from the total stellar mass of 
MD41, assuming half of it is enclosed within a radius of R e , 
is log(£(M*, r < R e )) = 7.6; for this fairly inclined system, this 
increases significantly when considering the projected surface 
density to log(S(M*,r < R e )) = 8.1, virtually the same as in- 

17 The two blue sources offset from the main part of MD 41 to the north- 
west and southeast that are detected in the ACS image but not in the NIC2 
image are included in this large aperture, and so contribute to this estimated 
integrated color. 



ferred on average at individual locations across the galaxy. 
Compared t o the propert i es der ived for massive 1 .5 < z < 2.5 
galaxies by iFranx et al. I (120081 see also Figure [TTTi. MD41 
lies at the low surface mass density end, consistent with its 
high specific SFR of 24 Gyr -1 and blue estimated integrated 
rest-frame ({/-B)ab ~ 0.5 color. 

The total stellar mass obtained by summing up the inferred 
stellar mass of individual pixels with S/N > 3 in the i'sh and 
Z/160 band images is 7.3 x 10 9 M©, or 95% of the value de- 
riv ed fr om the modeling of the integrated optical/near-IR SED 
(§ [33). The sum of the //160 band flux densities over the 
same pixels corresponds to 65% of the total as measured in 
a 3 "-diameter aperture (§ HTZT i. Multiplying the S/N > 3 
pixels' luminosity densities by the M*/L^ st from the best- 
fit SED model of MD41 results in a stellar mass of only 
5.0 x 10 9 M Q . This difference with respect to the spatially- 
resolved estimate reflects the difference between the pixel- 
averaged M*/U^ x and that from the integrated SED noted 
above. It suggests the total SED-derived stellar mass may 
be underestimated by sa 30% because dustier and/or more 
evolved, redder stellar populations are under-represented in 
the effectively luminosity-weighted integrated rest-UV and 
optical photometry of MD 41 . A similar effect between stellar 
masses computed by integrating spatially resolved mass maps 
and by mo deling unresolved photometry was found for nearby 
galaxies (Zibett i. Chariot. & Rixl 120 09). This effect may be 
more important for MD 41 (possibly also for BX 482) because 
the SED fits rely more heavily on the rest-UV photometry 
than for the other galaxies with better rest-optical coverage 
(three near-IR bands compared to two for MD 41 and one for 
BX482). 

Our analysis of ACS ;'si4 and NIC2 band imaging 
of MD41 shows modest variations in M*/L" st ratio, within 
a factor of ~ 1.5 over most of the galaxy and on spatially- 
resolved scales as small as pa 1 .2 kpc. In addition, the results 
indicate that the 1 .6 fim (rest-frame pa 5000 A) is a fairly good 
tracer of the stellar mass distribution of MD 41 . The spatially- 
resolved color-based mass map suggests nevertheless that 
luminosity-weighting biases may have led to a pa 30% under- 
estimate of the total stellar mass based on modeling of the 
integrated SED of MD 41. It will undoubtedly be interesting 
to investigate this effect more systematically with larger high 
z galaxy samples from combined high-resolution optical and 
near-IR imaging in future studies. 

7.2. The Sizes of our SINS NIC2 Sample Galaxies 

Size is a fundamental property of galaxies. The redshift 
evolution of disk sizes and scaling relations involving size 
(e.g., size-velocity, size-mass) set sensitive constraints on 
the coupling between the baryonic and dark matter, on 
the mechanisms through which disks acquire their angular 
momentum, a nd on the physical processes that drive galaxy 
growth (e.g.. IWhite & Reesl \Wm iFall & Efstathioul1l980l: 
Blumenthal et al. 1984; Dalcanton, Spergel, & Summers 
1997J IMo. Mao. & White! 11991 iNavarro & Steinmetz 
2000). Size determinations, however, can be affected by 
various factors, and reconciling observed siz es with theo- 
retical predictions remains a challeng e (e.g.. IDutton et al. I 
| 2007t IDutton & van den Boschl l2009bl IDutton et all 1201 lb 
Some rvilleetal.l 120081: ISaleset al.l 120091: iBurkert et all 
120101) . In the following, we consider whether the sizes 
derived for our NIC2 sample are affected by their particular 
morphologies and interpret them in the broader context of 
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the size distribution of the general z ~ 2 galaxy population. 
Our sample has high quality data that allow, for the first time, 
a reliable and consistent comparison of rest-frame optical 
continuum and Ha sizes, and we discuss implications for 
disk formation at high redshift. 

7.2.1. Rest-frame Optical Continuum and Ha Sizes 

It is clear from our NIC2 images that galaxies at z ~ 
2 can have complex, clumpy, and irregula r morphologies 
even i n the rest-frame optical (see also, e.g.. iPapovich et aD 
20051: iToft et al. Il2007h IZirm et al. Il2007l: iDasyra et al.ll2Q08h 
Kriek et al71l2009l) . Depending on the importance of the sub- 
structure, the overall radial light profiles may or may not be 
well reproduced by a simple smooth model (see Figures HI 
[5] O, and this substructure may also contribute to the ob- 
served scatter in light concentration index. This can compli- 
cate the measurement and interpretation of sizes in a broader 
context. For our NIC2 targets, we estimate in Paper II that the 
clumps identified in rest-frame optical continuum light con- 
tribute individually typically only ~ 2% of the total emis- 
sion of the host galaxy; their total contributions range be- 
tween ~ 10% and ~ 25% (15% on average), depending on 
the object. Therefore, for these six galaxies, the impact of 
the substructure associated with the clumps, spread across the 
galaxies, probably does not affect the size measurements in 
an important way. 

As we found in § 14.1.21 the kinematically identified disks 
have global radial light profiles that deviate appreciably from 
pure exponential disks, with best-fit Sersic indices n sa 0.15- 
0.6, or n sa 2 for the AGN source BX663. The impact of 
fitting exponential disks to our NIC2 targets is, however, rela- 
tively modest. Excluding the major merger BX528, single- 
component Sersic model fits to the 0."145 resolution H\a) 
band images with fixed n = 1 yield the following changes: 
the effective radii mostly increase, with mean and median 
R e (n = l)/i?e(« free) « 1.15, and range ofO.8-1.45. Alln < 1 
sources have R e (n = 1) > R e (n free) (with mean/median factor 
wl.2 and range of 1 . 1 - 1 .45), whereas BX 663 with n = 2 has 
R e (n = 1) = 0.8 x R e (n free). For fits to the Ha maps at sa 0"5 
resolution, where all five disks have a best-fit n < 1 (Table|5]l, 
the R e (n =1) are larger compared to R e (n free) by a factor of 
1.1 on average and medi an, ranging from 1.0 to 1.14. 

As discussed in § 16.31 the rest-frame optical sizes derived 
for our massive and actively star-forming targets, though 
large, are fully consistent with the trend of increasing ef- 
fective radius at higher specific SFR and lower stellar mass 
surface density observed among the general massive galaxy 
population at high redshift. While the objects studied in this 
paper include some of the larger galaxies from the SINS sur- 
vey (Figure [3), they lie well within the distribution of sizes 
derived for 1.5 < z < 2.5 galaxies at very similar rest-frame 
optical wavelengths (close to g-band) and in the relevant stel- 
lar ma ss range of log(M* [M©]) = 9.9 - 1 1 .0 by iFranx et aU 
(2008), which span Z? e ,tirc ~ 0.5 - 8 kpc. Our targets have 
^e.circ = 2.6—4.3 kpc, which places them within the upper half 
of the full distribution. The substantial scatter by roughly an 
order of magnitude in rest-frame optical (circularized) effec- 
tive radius at fixed stellar mass is seen at all redshifts up to 
z ~ 2.5 and is attributed, at least in part, to the wide range 
of star formation rates and histories among galaxies (e.g. , 
Shen et al. 120031: [Kauffmann et al.l2003 T: lTrujillo et aD l2006l: 
Franx et al. 112008b IToft et al. 1120091 and see Figure [TBI This 
scatter as well as the relationships between size, star forma- 
tion activity, stellar mass, and stellar mass surface density 



mean that sample properties and selection biases need to be 
considered for consistent comparisons with other galaxy sam- 
ples as well as theoretical predictions. 

The rest-frame sa 5000 A continuum probed by our NIC2 
data is more sensitive to lower-mass, longer-lived stars mak- 
ing up the bulk of the stellar mass while the Ha emission is 
dominated by the young, short-lived stars. One might expect 
to be able to discern with these data possible differences be- 
tween the distributions of stellar mass and of star formation, 
if any. For all our SINS NIC2 galaxies, however, we find 
very close agreement between the parameters derived from 
the 7/i6o band and Ha maps, with on average the same ef- 
fective radius and Sersic index to within 3% and 0.08, re- 
spectively (§ 15.21 and Figure [T3l. This similarity in size and 
global profile shape is probably not due to beam-smearing ef- 
fects, which we found to be not very significant for our targets 
based on the Z/ieo band images at the original NIC2 resolution 
and convolved to the typical sa 0."5 resolution of the SIN- 
FONI data as well as on the seeing-limited and AO-assisted 
Ha maps of BX 482. 

These results, albeit for a small sample, are interesting to 
consider in the context of models of disk formation, where a 
signature of "inside-out" scenarios for the growth of stellar 
disks is that the half-mass radius of the recently formed stars, 
"/? e (SF)," is larger than the half-mass radius from the total 
stellar ma ss, "RejM*)." In par ticular, in the semi-ana lytical 
models of Dutt on et al.l ( 1201 U see also lDuttonll2009al) . the 
R e (SF) is on average a factor of sa 2 larger than R e (M+) 
for z ~ 2 disks. Us i ng cos mological hydrodynamical sim- 
ulations, ISales et aTl d2009h find a similar difference by a 
factor of « 1.8 at z~ 2. To account for rad ial variations 
in mas s-to-light ratio across their model disks, iDutton et al.l 
(1201 11) also compute half-light radii in several rest-frame op- 
tical and near-IR bands, and find that the differences be- 
tween half-light radii and half-mass radius of recent star for- 
mation decreases towards bluer bands. For the V band, the 
bandpass that most closely matches the rest-frame emission 
pr obed by the Z /i6p observations of our targets, the models 
of IDutton etaTI imply on average R e (SF)/R e (V) w 1.26. In 
these models, "recent star formation" refers to stars formed in 
the last sa 30- 100 Myr, and we assume here that it is traced 
reliably by Ha. Excluding the major merger BX 528, the rest- 
frame optical radii derived from the NIC2 images matched to 
the PSF of the SINFONI Ha data in Table |5] give on average 
R e (Ha) / R e (Hi(,o) = 1.06 (and 1.02 for the median), ranging 
from sa 0.95 (MD41, BX482) to sa 1.15 (BX663, BX610). 
These results for fits where n was left as a free parameter 
are appropriate for comparison with the IDutton et al. I (1201 ll) 
models, in which disks are not formally exponential. For com- 
pleteness, the same comparison based on fits with n fixed to 
1 gives an average R e (Ha) / R e (H\6o) ratio of 1.10 (median of 
1.08), in the range sa 1.0-1.2. Formally, the uncertainties of 
individual objects would allow ratios of sa 1.25 but the trend 
of ratios sa 1 .05 - 1 . 1 (with standard deviation of 0. 1) for the 
ensemble of our five z ~ 2 disks indicates little, if any, system- 
atic difference between the Ha and rest-frame optical sizes. 

We conclude from the considerations above that the size 
measurements of our SINS NIC2 targets appear robust against 
possible biases due to substructure, or to specific assumptions 
on the light profile as tested with fixed n = 1 fits to our disks. 
We estimated the impact of substructure or of adopting n = 1 
models to all be within typically ~ 15%. Based on the NIC2 
and SINFONI data of the five disks among our sample, we 
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do not detect a clear effect whereby the distribution of the ac- 
tively star-forming regions would be characterized by a larger 
effective radius relative to that of stars tracing more closely 
the bulk of the stellar mass. It would be premature to gener- 
alize this finding from a small number of galaxies and draw 
conclusions on the "inside-out" scenario of disk formation, 
but our high quality NIC2 and SINFONI data provides a first 
set of direct empirical constraints at z ~ 2 for theoretical mod- 
els. The sizes of all our six galaxies, with R e pa 4-6 kpc (or 
^e.circ ~ 2.5 -4.5 kpc) are not surprisingly large when placed 
in the broader context of the massive z ~ 2 galaxy popula- 
tion, taking into account their properties in conjunction with 
the existing trends of larger sizes at high specific SFR and low 
stellar mass surface density. 

7.2.2. Possible Caveats 

There are several possible caveats to the analysis presented 
in the previous subsection. In particular, our simple Sersic 
model fits do not account for multi-component structure. As 
a specific example, if there is a significant bulge, the R e from 
single-component model fits will tend to underestimate the 
R e for the disk itself. The Hk,o band radial light profiles 
do not indicate the presence of a significant bulge compo- 
nent in our disks (Figure |6j, and this is further supported by 
the low Sersic indices and relatively small Gini coefficients 
(Figures [T3l and [T4l). The exception may be BX663 but it 
is unclear whether, and by how much, the AGN in this ob- 
ject may contribute to the central emissi on pea k. The two- 
component fits for BX 663 explored in § 14.1.21 imply a very 
minor fraction of 6% of the total Hk,o band light for a pu- 
tative nuclear point source but a more important 35% for a 
putative bulge. In both cases, however, the R e of the disk 
component is not significantly different from that of the best- 
fit single-component model (with n = 2). Also, as we deter- 
mine in Paper II, the clump closest to the center in each of the 
other four disks, which could be interpreted as a young bulge, 
makes a marginal contribution of < 5% to the total rest-frame 
pa 5000 A light, and even only pa 13% for the more prominent 
central peak of BX 663(3 It thus appears unlikely that the 
sizes inferred for the disks in our SINS NIC2 sample would 
be significantly affected by the presence of a bulge, neglected 
in our Sersic model fits. 

Another caveat is that the effects of dust obscuration are ig- 
nored. Extinction may significantly modify the surface bright- 
ness distributions, and hence the best-fit sizes. In particular, 
higher extinction towards the central parts of the sources will 
lead to shallower observed inner profiles. For BX610, our 
seeing-limited SINFONI H and K band data are of sufficient 
S/N to constrain the radial variations of the H /3/Ha flux ra- 
tio (IP. Busc hkamp et" al71l201 1. in pr eparation). The Balmer 
decrement suggests indeed an increase in visual extinction 
by pa 1 mag in the inner pa 5 kpc, so that the intrinsic inner 
light profile could be steeper and the effective radius could be 
smaller. Moreover, differential extinction towards the H II re- 
gions and the stars could affect the relative rest-frame optical 

18 From detailed modeling of the Ha kinematics of SINS disks with high 
quality SINFONI data, including four of our targets, Genzel et al.] <2008D 
found a trend of increasing central dynamical mass concentration with stel- 
lar evolutionary stage, which could be interpreted as evidence for growing 
bulges/inner disks. The inferred ratio of enclosed dynamical mass within 
3 kpc and 10 kpc ranges from < 0.15 for MD41, to 0.2 for BX482, and 0.4 
for BX389 and BX610. Disentangling the relative contributions of a bulge 
and the underlying disk from these data is uncertain but these estimates would 
provide conservative upper limits to bulge-to-disk mass ratios. 



and Ha effective radii. The existing data for our NIC2 targets 
do not allow us, however, to constrain on ~ kpc scales the 
effects of spatial variations in dust extinction and of differen- 
tial obscuration between the H II regions and the bulk of the 
stellar populations. 

Finally, the depth of the Hiw band and Ha maps may also 
influence the best-fit Sersic model parameters, including the 
n index and R e . Our SINFONI and NIC2 data sets are fairly 
deep in term s of wha t can be achieved with these instruments 
(see §§ 12.21 and I2.6K and the radial light profiles are traced 
to > 2 x R e (Figures [6] and [T2l . Nonetheless, the 3 a sur- 
face brightness limits correspond to fairly high Ha and rest- 
frame g band luminosities of Lu a pa 4.3 x 10 39 ergs" 1 kpc" 2 
and L™ 1 pa 6.7 x 10 8 L g , Q kpc" 2 at the median z = 2.2 of our 
NIC2 targets, and it is conceivable that fainter extended emis- 
sion in the outer parts is undetected in our data. 

7.3. Relation Between Kinematics and Morphologies 

Arguably, one of the key aspects in the study of high red- 
shift galaxies is a characterization of the mass assembly pro- 
cess, and, in particular, pinpointing evidence for the occur- 
rence of major mergers and interactions or, in their absence, 
for the establishment of disk-like systems, which are a pre- 
diction of models of galaxy formation. There are several 
inherent difficulties in studies based on morphologies, in- 
cluding bandshifting, intrinsically complex structure, and un- 
certainties in the timescal es over which key merger signa- 
tures are observable (e.g.. iToft et al.l|2007t [Elmegree n et al.l 
120071 120091; lLotz et al. ll2008l l2010allbh . Therefore, kinemat- 
ics appear to be a prerequisite for a robust distinction al- 
though there are possible complications as well (for instance, 
particular m erger co nfigurations mimi c king disk-like kine- 
matics; e.g.. lLaw. Steidel. & Erbll2006t iShapiro et al 1 1 20081: 
[ Robertson & Bullocki 120081; IForster Schreiber et al.l 120091 
Epinat et al. 1 120 1 Of) . Morphologies and kinematics together 
make for a powerful combination, mitigating the respective 
potential limitations as well as providing a more co mplete pic- 
ture of individual galaxies (see also lNeichel et al.ll2008l) . 

The kinematic classification of our six SINS targets, all 
disks except BX528, was first based on SINFONI data 
sets taken with the 0."125 pixel" 1 scale, at pa 0."5 or pa 
4.1 kpc resolution (0."39 and 3.2 kpc for the NGS-AO 
data of BX663). Initially based on qualitative assess- 
ment of b oth the Ha velocity field and velocity dispersion 
map (IForster Schreiber et al.ll2006bl) . the classification into 
disk- and merger-like systems was qu antitatively confirme d 
through the application of kinemetry ( Sha piro et al.l [2008), 
with the exception of BX663 due to S/N limitations out- 
side of the inner regions. To date, BX482 is the only tar- 
get among our NIC2 sample followed-up with AO and the 
0."05 pixel" 1 scale, with effective resolution of 0"17. Its kine- 
matic classification as a disk is further confirmed with the 
higher resolution data, although more detail in the Ha emis- 
sion and kinematic maps is resolved. Fourteen other galax- 
ies from the full SINS survey were classified by kineme- 
try, of which four have follow-up AO obse r vations at ~ 
0"15 - 0."20 resolution (IShapiro et a ll 12001 IGenzel et al.l 
[2001 IForster Schreiber et al.l 120091; IGenzel et al.l 120101; 
P. Buschkamp et al. 201 1, in preparation). In all cases, as for 
BX 482, the kinematic classification is the same whether it is 
based on the seeing-limited or the AO data sets. 

At first glance, the NIC2 //160 images of our six targets 
at 0."145 resolution (pa 1.2 kpc) reveal rich sub-structure at 
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all surface brightness levels. Qualitatively, multiple bright 
compact components can be interpreted as distinct interact- 
ing/merging units, asymmetries in the more diffuse emission 
or distorted i sophotes on large scales can be attribute d to tidal 
features (e.g. lLotz et"aLl l2004; Con selice et al.ll2005l) . On the 
other hand, these characteristics may also result from inter- 
nal dynamical processes within gas-rich disks, such as bar- 
and spiral-like features from non-axisymmetric perturbations 
or luminous kpc-sized clumps resulting from disk instabilities 
and tidal-like features induce d by the mutual interactions be- 
tween massive clumps (e.g.. iBournaud et al. 1120071 120081) . It 
may therefore not be immediately obvious whether a system 
is a disk or a merger f rom si mple visual inspection. Quan- 
titatively, as found in § 14.1.21 all the kinematically classified 
disks among our NIC2 sample have Sersic indices n < 2 (and 
even < 1 when excluding the AGN source BX 663), consis- 
tent with structurally late-type, disky systems. 

Among our NIC2 sample, BX610 appears to be a par- 
ticularly striking example where the rest-frame optical mor- 
phology reflects tell-tale signatures of internal d ynam ical pro- 
cesses. Indeed, its appearance, as outlined in § 13.11 exhibits 
evidence for a bar and features reminiscent of spiral arms, and 
there is no evidence for any nearby companion in our NIC2 
and SINFONI data, nor of recent merging. This galaxy is 
one of the most robust exam ples of rotating disk kinemat- 
ics among the SINS survey dForster Schreiber etal.l[2 006b; 
IGenzel et all 120081: IShapiro et al.ll2008t ICresci et al.(l2009l) . 
with the Ha velocity field exhibiting the classical "spider di- 
agram" of isovelocity contours, a nearly axisymmetric and 
centrally peaked velocity dispersion map, and a close corre- 
spondence of the kinematic center and major axis with the 
geometrical center and morphological major axis (from both 
the rest-frame optical and Ha emission). BX610 shows sev- 
eral prominent clumps and direct observations of CO molec- 
ular gas indicate a very high gas to baryonic mass fraction of 
^55% dTacconi et al. I2.010T) . All these properties suggest that 
BX610 may be one such apparently isolated galaxy where 
clumps formed from fragmentation in a Toomre-unstable 
gas-rich disk and where effic ient secular proces ses drive 
its early dynamica l evolution (|Genzel et al. 1 120081 see also 
Carollo et al 1120071 ^IBournaud et al ll2007tlDekel et al. I l2009t 
Agert z et al. Il2009t iCeverino et al.ll2010lh While BX 610 is a 
particularly remarkable case, the other disks among our NIC2 
targets as w ell as other high z samples studied at high res- 
olution (e.g.. IBournaud et aL]l2008t lElmegreen & Elmegreen 
120051: lElmegre enet al Il2009t IGenzel et al. 1 120 1 PL and Paper 
II) show prominent clumps and properties consistent with the 
above theoretical scenario. 

While the majority of objects in our SINS NIC2 sample 
are classified as disks, it is noteworthy to highlight examples 
of how the unique combination of Ha kinematics and rest- 
frame optical morphologies sheds light on galaxy interactions. 
For BX528, the Ha kinematics clearly indicate its merger 
nature even at the ~ 0"6 seeing-limited resolution of the 
SINF ONI data dForster Schreiber et al. Il2006bl: IShapiro et al.l 
2008). The reversal in velocity gradient across the system 
could be produced by a counter-rotating binary merger. In the 
NIC2 image, the system is resolved into two main "SE" and 
"NW" components separated by about 1" (8 kpc) in projec- 
tion, and connected by a fainter "bridge" of emission. Single- 
component Sersic model fits to the //j6o map, either at the 
original NIC2 resolution or convolved to the lower resolution 
of the SINFONI observations, imply a shallow inner light pro- 
file with n ~ 0.2, evidently caused by the lower surface bright- 



ness between the SE and NW components. As seen in § 13.11 
at 0."6 resolution, both the Ha and smoothed //i6o band mor- 
phologies appear rather regular and the SE and NW compo- 
nents are not discerned, making it impossible to distinguish 
from disk-like systems. Moreover, BX 528 does not stand out 
compared to the disks in our sample in the G-'J — Moo pa- 
rameter space based on the high resolution NIC2 image (Fig- 
ure 0. Thus, for BX528, both the Ha kinematics and the 
high resolution rest-frame optical imaging are necessary to 
fully characterize the nature of this system. 

Even if not as significantly as BX 528, BX 389 also presents 
evidence for galaxy interaction (Figure ITOb . The Ha kine- 
matics at 0"5 resolution show fairly regular motions consis- 
tent with rotation in a highly inclined disk. The various line 
maps from the K and H band seeing-limited SINFONI data 
all exhibit an extension to the south of the center, and the 
rest-frame optical line ratios ([N II]/Ha, [O III]/H/3) are con- 
sistent with either a low-metallicity companio n or shock ex- 
citation in a starburst-driven galactic outflow ( Lehnert et al. I 
120091: IP. Bu schkamp et al7l l201 l.Tn preparatipn|). This exten- 
sion does not appear kinematically distinct from the main part 
of the source, with line-of-sight velocities similar to those 
along the major axis of the main body of the galaxy, on the 
same side and at the same projected radii. The NIC2 H\a) 
band imaging reveals that this is a separate component, at 
a projected distance of 5 kpc below the galactic plane of 
BX389, supporting the nearby satellite interpretation. For 
BX389, therefore, the high resolution near-IR imaging re- 
solves the companion from the main part of the source and 
the kinematics confirm the physical association between the 
two components. 

Assuming that the observed 1 .6 /im lig ht trac es closely the 
stellar mass in our NIC2 targets (see § 17.1.21 1. we can set 
constraints on the mass ratios of the different components 
in BX528 and BX389. Variou s est imates based on the two- 
component Sersic model fits (§ 14. j} , on simple aperture pho- 
tometry, and on the fractional light contributions of the associ- 
ated "clumps" (presented in Paper II) imply a mass ratio in the 
range 1.2:1 to 1.9:1 for BX528-SE relative to BX528-NW, 
and 6:1 to 20:1 for the main part of BX389 relative to the 
southern companion. From the stellar mass ratio, BX 528 is 
a nearly equal-mass major merger, in agreement with the sig- 
natures and classification from the Ha kinematics. In con- 
trast, the southern companion of BX 389 represents a minor 
merger, unlikely to have strong impact on the large-scale kine- 
matics of the primary component and consistent with the lack 
of substantial disturbances in the kinematics expected for ma- 
jor mergers. 

We conclude from the above that the Ha kinematics and 
rest-frame optical continuum morphologies of our NIC2 sam- 
ple draw a consistent picture for the nature of each target and 
are very complementary. The kinematics, even at ~ 3 - 5 kpc 
resolution, are sufficient (and in some cases necessary, as 
for BX528) to reliably distinguish between disks and ma- 
jor mergers. We caution that there will be specific cases 
where a system's configuration will lead to ambiguous or 
erroneous classificat i on, also dependi ng on the data quality 
dShapiroet al.l 120081: lLawet all 120091) . but this does not ap- 
pear to be the case for the objects studied in this paper. Mor- 
phologies (as well as kinematics) at higher ~ 1 kpc resolu- 
tion then reveal further details and are important to identify 
minor mergers, and more subtle or s maller-scale signatures 
of internal dynamical pr ocesses (e.g. IBournaud et aU 120081; 
IGenzel et al. 1120081 l2010t) . 
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8. SUMMARY 

We have presented //5T/NIC2 F160W imaging for a 
unique sample of six z ~ 2 star-forming galaxies with 
VLT/SINFONI integral field spectroscopic measurements. 
Our analysis has been based on the combined NIC2 and SIN- 
FONI seeing-limited and AO-assisted datasets, along with 
multi-wavelength ground-based photometry and, in one case, 
HST/ACS imaging. Using NIC2 imaging alone we have car- 
ried out galaxy structural analyses, using both parametric and 
non-parametric approaches, and placed the structural prop- 
erties of our SINS NIC2 sample in the larger context of the 
correlations between morphology and star-formation activ- 
ity present at z ~ 2. However, it is the full combination of 
datasets presented here that has enabled investigations of un- 
precedented detail into the dynamical nature of our targets, 
and the importance (or lack thereof) of mergers of different 
mass ratios; the relative galaxy sizes traced by recent star for- 
mation and stellar mass; and spatial variations in the stellar 
M/L and stellar mass surface density. 

The main scientific conclusions of this paper can be sum- 
marized as follows. 

1 . The rest-frame optical continuum morphology of our tar- 
gets yields classifications consistent with those based on ion- 
ized gas kinematics. Those systems with disk-like kinemat- 
ics (the majority of our sample) are characterized by disk-like 
structures in rest-frame optical light. Furthermore, for the one 
target (BX528) with kinematics indicative of binary counter- 
rotating elements, the NIC2 image is clearly resolved into 
two, discrete components separated by w 1", or 8 kpc. While 
BX 528 may provide evidence for a major merger, the mor- 
phology and kinematics of BX 389 provide evidence for a less 
dramatic interaction. The SINFONI kinematics for BX 389 
indicate regular disk rotation, but both SINFONI Ha maps 
and NIC2 rest-frame optical imaging reveal a low-luminosity 
companion at roughly the same line-of-sight velocity, and a 
projected distance of 5 kpc — possibly a minor merger in ac- 
tion. 

2. We use both parametric (Sersic profile) and non- 
parametric {G,^,M2o) methods for characterizing the NIC2 
morphologies of the SINS NIC2 sample. The majority of 
these objects are characterized by shallow (« < 1) inner light 
profiles, with median effective radius of R & ~ 5 kpc. Signif- 
icant residuals are found with respect to the model fits for 
some objects due to the presence of multiple discrete clumps 
in the surface-brightness distributions, giving rise to an over- 
all irregular appearance. The G, Vf, and M20 values for the 
sample also reflect the clumpy appearance, with significantly 
higher M20 and ^ values than would be derived for axisym- 
metric, monolithic, and compact galaxies at similar redshift. 
The relatively low G values indicate the lack of any single, 
prominent bulge component. 

3. Given the small sample size, and the active levels of 
ongoing star formation in all of our targets, it is difficult to 
extract strong trends between morphology and other galaxy 
properties within the SINS NIC2 sample. However, a com- 
parison with other z ~ 2 samples selected using a rest-frame 
optical (observed Zf-band) or mid-IR (observed 24 ^m) limit 
reveals that our galaxies follow the broader trends of increas- 
ing size, M20, and ^P, and decreasing n, G, and stellar mass 
surface density, with increasing specific star formation rate. 
At the same time, the larger ensemble of galaxies reflects a 
considerable diversity in structure at fixed stellar mass. While 
these trends have been presented in other work using seeing- 



limited or lower-resolution NIC3 images, the sample with ad- 
equate resolution for detailed profile fits (i.e., NIC2) is still 
very small, and our objects represent an extension of the deep 
z ~ 2 NIC2 dataset to higher specific star formation rate. 

4. The data presented here allow for a comparison between 
rest-frame optical morphologies and those derived from trac- 
ers of active star formation such as Ha e mission and rest- 
frame UV light. So me theoretical models dSales et al. Il2009t 
iDutton et al. 1 1201 Tb predict that the size of the galaxy ob- 
tained from tracers dominated by young stellar populations, 
such as Ha emission, should be systematically larger than 
that traced by rest-frame optical light. For the SINS NIC2 
sample (excluding the one putative major merger, BX528), 
the Ha and NIC2 //160 band sizes reveal no significant dif- 
ferences. For the one object with AO-assisted SINFONI 
data obtained at comparable resolution to our NIC2 imag- 
ing (BX482), the non-parametric coefficients for Ha and 
rest-frame optical continuum emission are very similar. Fur- 
thermore, the average NIC2 G, and M20 coefficients for 
our UV-selected SINS NIC2 targets are very similar to those 
based on rest-UV ACS morph ologies for a large r sample of 
UV-selected galaxies at z ~ 2 (iPeter et al. 1 12007). MD41 is 
the only galaxy in our sample with both NIC2 F160W and 
ACS F814W imaging. The rest-UV morphology for this ob- 
ject is broadly similar to that in the rest-frame optical, al- 
though, in detail, its derived R e and n are slightly larger and 
lower, respectively. Overall, however, the parametric and non- 
parametric coeffients for this object are in agreement. For 
UV-selected z ~ 2 galaxies with active star formation, the 
rest-frame optical appearance thus seems very similar to that 
traced by rest-frame UV light and ionized gas. We have also 
established that the NIC2 band images are truly tracing 
continuum emission, and are not dominated by nebular line 
emission. 

The power of combining ^kpc-scale resolution rest-frame 
optical and UV imaging together with integral field spec- 
troscopy at seeing-limited and AO-assisted resolution has 
been demonstrated here for a small sample of galaxies. With 
such a dataset, it is possible to characterize simultaneously 
the past and present assembly of mass, and the history of star 
formation. A statistical sample of such measurements will 
provide crucial inputs into models of the formation and evo- 
lution of massive galaxies, and will be possible with current 
and planned facilities in space and on the ground. 
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APPENDIX 

A. DERIVATION OF THE RELATIONSHIP BETWEEN MASS-TO-LIGHT RATIO AND COLORS FOR MD 41 

In § 17.11 we combined the NIC2 Hiw band data with ACS i%u band imaging available for one of our targets, MD41, to 
investigate spatial variations in color, mass-to-light ratio, and stellar mass on resolved scales of 1.2 kpc. It is well known that 
obs erved colors and M/L ratio of stellar pop ulations are correlated, and this behavior can be used to infer stellar masses (see, 
e.g jRudnick et aU l2006: Zib etti et al. 1 12009, for applications to derive the cosmic stellar mass density or spatially-resolved mass 
maps of nearby galaxies). In this appendix, we present the derivation of the relationship we used to interpret quantitatively the 
hiA~H\aa color map of MD 41 . 

We employed iBruzual & Ch ariot (2003) models to generate the observed /si4-//i6o colors of stellar populations appropriate 
for z = 2.2, corresponding to the redshift of MD41. The colors were synthesized from the redshifted model spectra using the 
total transmission curves of the NIC2 F160W and ACS F814W bandpasses. We considered a variety of star formation histories 
(SFHs) ranging from single stellar populations ("SSP") with all stars formed instantaneously and no subsequent star formation, 
exponentially declining SFRs with e-folding timescales from 10 Myr to 1 Gyr, and constant star formation (CSF). Th e models 
were c omputed for a grid of ages from 1 Myr to 5 Gyr, and of extinctions with Ay from to 3 mag. We adopted the Chabrier 
(2003) IMF and the lCalzetti et al.l (120001) reddening law, and considered two different metallicities with Z = Zq and 0.2Z©. 

For every parameter combination (SFH, age, Ay, Z), the models provide the ratio of stellar mass to rest-frame optical luminosity. 
We consider here the "effective" M/L ratio, i.e., where the luminosity refers to the light that reaches the observer after being 
attenuated by dust. We further consider the luminosity at rest-frame wavelengths close to the range probed by our NIC2 F160W 
data of MD41. Theses choices are relevant because we ultimately combine the derived M/L map with the observed Hm,o 
band image (least affected by extinction and by light from young stars) to derive the stellar mass map. For our purposes, we 
approximated the rest-frame optical luminosity seen through the NIC2 F160W filter (mean wavelength of 1.61 /itn, range with 
peak-normalized transmission > 50% of 1.40- 1.80 /im, corresponding to rest-frame 5080 A and 4415-5675 A, respectively) 
with the model predictions through the SDSS g-band filter (mean wavelength of 4695 A and range of 4020-5330 A). The g- 
band is w 380 A bluer than the rest-frame range probed by our NIC2 data, but it is closer in mean wavelength and has larger 
overlap than the V band (with mean wavelength of 5510 A and range of 5010-5865 A). Over the entire range of SFHs and ages 
considered, the differences between M^/U^ 1 and M+/Ly st ratios are of at most a factor of two, and similarly over the range of 
extinction assumed. Metallicity has a very minor impact in comparison. These differences in M/L ratios constitute however 
extremes as they are for the full range of model parameters. More realistically, the fairly narrow range in observed !8i4-fli60 
colors across MD41 suggest only modest variations in age and extinction. With the assumption of ages ~ 10 7 - 10 s _yr and 
Ay ~ 1 mag, more consistent with the results from the modeling of the integrated optical to near-IR SED of MD 41 in § 13.31 the 
M+/L™ 1 ratio is « 20% lower than the M*/Ly st ratio, with little dependence on the assumed SFH. Since our NIC2 F160W data 
probe a rest-frame wavelength range intermediate between the g and V bands, the numbers above represent upper limits on the 
offset between M+/L™ st we use and the actual stellar mass to rest-frame optical luminosity ratio from our data. 

Figure [l9]plots the observed igi4 — #160 colors versus M+/L^ sl ratios derived from the models for solar metallicity. Curves with 
different colors correspond to different SFHs. Stellar age increases along the curves from blue to red isu—Hieo colors, and from 
low to high M/L ratios. Different line styles are used to plot models with different extinction (shown for Ay =0, 1,2, and 3 mag 
in the Figure). The main effect of increasing extinction is to shift the model curves along a direction that is roughly parallel to 
the age tracks. The model curves occupy a fairly well defined locus in the observed /si4-#i60 versus M+/L™ st parameter space, 
reflecting the strong degeneracy between stellar age, extinction, and SFH in these properties. This degeneracy can be exploited 
to derive the mean effective M/L ratio for a given observed color, which, multiplied by the observed luminosity, yields a stellar 
mass estimate without requiring any knowledge about the actual SFH, age, or extinction. 

We computed the "mean" relationship by taking the average of the minimum and maximum values in log(M*/L 1 f est ) over all 
model curves within equally-spaced bins of width 0.1 mag in i%u— Hi^o colors. By proceeding in this way, we reduced the 
impact of the discreteness of our grids of SFHs and extinction. Over the ranges in colors where models for all SFHs and all 
ages for a fixed extinction value are represented, the average of the minimum and maximum log(M + /L' ; est ) is very close to the 
weighted mean over all models, where the weights are obtained from the time span of each model within a given color bin; the 
differences are at most 0.08 dex in log(Mi/L" st ). Varying the Ay value has little impact since extinction effects are mostly along 
the relationship. The mean relationship for solar metallicity is plotted in Figure [19] (white-filled circles and thick solid black 
line), with error bars corresponding to the standard deviation of all models within each color bin. The relationship derived for 
1 /5 solar metallicity is also shown for comparison (grey-filled circles and thin solid black line), but is very similar to the one 
for solar metallicity. For our analysis of MD41, we adopted the /8i4-#i60 versus log(M*/L" st ) relationship derived for Z = Z Q 
models. The standard deviation of the models suggest that the relationship is accurate to 0.15 dex in log(M*/L r ? est ) over the range 
of observed colors across MD 41 (indicated by the shaded histogram in Figure [19l>. 
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Kpc-scale rest-optical imaging of z ~ 2 galaxies 



TABLE 1 

Galaxies observed 



Object 




Kinematic Type b 




G 


n 


J 


H 160 










(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Q1623-BX528 


2.2683 


Major merger 


24.52 ±0.1 8 


23.81 ±0.14 


23.56±0.13 


22.44 ±0.17 


22.33 ±0.06 


21.61 ±0.20 


Q1623-BX663 C 


2.4332 


Disk 


25.40 ±0.24 


24.38 ±0.17 


24.14±0.15 


23.41 ±0.30 


22.79 ±0.10 


21.78±0.21 


SSA22a-MD41 


2.1704 


Disk 


24.81 ±0.11 


23.50 ±0.06 


23.31 ±0.05 




22.64 ±0.05 


22.30 ±0.36 


Q2343-BX389 


2.1733 


Disk 


26.39 ±0.39 


25.13±0.19 


24.85 ±0.16 


23.82 ±0.24 


23.11±0.10 


22.04 ±0.25 


Q2343-BX610 C 


2.2103 


Disk 


24.67 ±0.21 


23.92 ±0.12 


23.58±0.11 


22.35±0.13 


22.09 ±0.06 


21.07±0.13 


Q2346-BX482 


2.2571 


Disk 


24.44 ±0.1 2 


23.54±0.10 


23.32 ±0.09 




22.34 ±0.07 





NOTE. — Magnitudes are given in the AB system. Optical, J and A" s magnitudes are taken from Erb et al. 2006 (with the 
exception of SSA22a-MD41 for which the K s magnitude was measured on NTT/SOFI imaging available from the ESO archive). 
The Hi6o band magnitudes were measured from the new NICMOS/NIC2 data presented in this paper. 

a Systemic vacuum redshift derived from the integrated Ha line emission in the SLNFONI data {Rj rster Schreiber et al. 2009). 
b Classification according to the Ha kinematics from SINFONI ( Shapiro et al. 2008; Forster Schreiber et al. 2006b, 2009). 
c Spectral signatures of an AGN are detected in the optical and near-IR spectrum of 01623 -BX663 and the large-scale Ha 
kinematics suggest the host galaxy is a low-inclination rotating disk (Erb et al. 2006; Forster Schreiber et al. 2006b. 2009). The 
Spitzer/MIPS fluxes of Q1623-BX663 and Q2343-BX610 are indicative of obscured AGN activity JReddv et al.H2010h . 



TABLE 2 

LOG OF NICMOS/NIC2 OBSERVATIONS 



Object 


Date observed 


P.A. a 


3 a Limit in d = 0.22" b 






(degrees) 


(mag) 


Q1623-BX528 


2007 Apr 27 


-5.43 


28.07 


Q1623-BX663 


2007 Apr 28-29 


-20.43 


28.17 


SSA22a-MD41 


2007 May 07-08 


30.13 


27.98 


Q2343-BX389 


2007 Jul 10 


16.91 


28.08 


Q2343-BX610 


2007 Sep 24 


-45.43 


28.14 


Q2346-BX482 


2007 Jul 18 


4.57 


28.22 



NOTE. — The total on-source integration time is 10,240 s for all targets. 
:l Position angle of y-axis (in degrees East of North). 

b Effective 3 a background noise limiting magnitudes (AB system) in a "point-source" circular aperture with 
diameter d = 0''22, or 1 .5 x FWHM of the NIC2 PSF as measured from stellar profiles in the data. The limiting 
depths were determined from the analysis of the noise properties of the reduced data presented in § 12.61 



TABLE 3 

Stellar and dust properties of the galaxies 



Object 


Age 


Av 






SFR 


sSFR b 




(Myr) 


(mag) 


(10 10 M Q ) 


(mag) 


(Moyr" 1 ) 


(Gyr 1 ) 


Q1623-BX528 
Q1623-BX663 
SSA22a-MD41 
Q2343-BX389 
Q2343-BX610 


2750+ % 
3 -21,10 

2500+ 147 

50+ 31 
2750+224 


0.6±0.2 
0.8 ±0.2 
1.2±0.2 
1.0±0.2 
0.8 ±0.2 


6.95+ 017 
o.fu_ 2 28 
77+°- 10 
4 j 2+0.77 

10 0^ 


-22.90+™ 4 , 
-22.68+SSl 

99 Q7+0.04 

_21. 94™ 09 
-23.10+^ 


42+29 
49+ 1 3* 

185 -70 
25+' 7 
60+f 
80+« 


6+ 1 7 
07** 

u '-0.2 

24+' 

6+ 1 2 
u -°-o.i 

0.6±0.2 


Q2346-BX482 


Q9 1 +485 

3* l -M 


0.8 ±0.2 


1 84+°™ 


-22 65+^ : ^ 
ZZ - OJ -0.08 


4 3+3.0 



NOTE. — Results from SED modeling of the optical and near-IR photometry, with best-fitting star formation 
history having a constant SFR for all targets. The formal (random) fitting uncertainties are g iven, derived from 
the 6 8% confidence intervals based on 200 Monte Carlo simula tions for the default set of iBruzual & Charlotl 
2003 models with solar metallicity, a Chabrier 2003 IMF, and the Calzetti et al. 2000 reddening law; systematic 
uncertainties (from SED modeling assumptions) are estimated to be typically a factor of ~ 1.5 for the stellar 
masses, ±0.3 mag for the extinctions, and factors of ~ 2-3 for the ages as well as for the absolute and specific 
star formation rates. 

a Rest-frame absolute V-band magnitude, uncorrected for extinction. 

b Specific star formation rate, i.e. the ratio of star formation rate over stellar mass. 
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TABLE 4 

Best-fit Structural Parameters from the NIC2 H l60 Band Images 



Object 


Re 

(kpc) 


11 


b/a 


P.A. 

(degrees) 




Single-component models 


Q1623-BX528 
Q1623-BX663 
SSA22a-MD41 
Q2343-BX389 
Q2343-BX610 
Q2346-BX482 


t.Jt_ 079 

5 69+ 20 

-Q-ia 

c Q3+0.I7 

4.44 ±0.08 

6 22+ 013 


0.16±0.01 

2 00+ 211 

54+ 007 
-8 0S 

36 + ° i " 
u - -0.05 

0.57 ±0.04 
0.14±0.01 


0.31 ±0.01 
73+ 004 

"•' J -0.02 

0.31 ±0.01 
0.30 ±0.01 
0.56 ±0.01 
0.48 ±0.01 


-27.8 + [' 2 

-21.8^9 
36.5±0.2 

-48.6+° 3 
-60.8!}] 


Two-component models 


Q1623-BX528 SE 
Q1623-BX528NW 
Q2343-BX389 Main 
Q2343-BX389 S 


3 ,8+0.18 

3 =7+1. 98 
-Q-52 
5 RQ+0-08 

7 71+2.18 
z -"--0.69 


1 20+°' 15 

2 91+170 

n 07+o^f 
u - JZ -0.04 

3 03+ 2 64 


35-° m 

J -0.03 

051 +o.06 

U - J1 -0.05 

0.27 ±0.01 
20+ 35 


^K).5 +3 - 4 
-10.9^| 
-49 n + ° 5 
-41 5+ 22 



NOTE. — Structural parameters derived from Sersic profile fits to the two-dimensional rest-frame ~ 5000 A surface 
brightness distribution from the NIC2 Hk,o band images. Single-component models were fitted to all targets. For BX528, 
two-component fits were also performed to account for the southeast ("SE") and northwest ("NW") components. Similarly 
for BX389, where the two components represent the main part of the galaxy and the small compa nion to the south, "Main" 
and "S ," respectively. The results are obtained from a series of 500 GALFIT runs as described in § 14.11 The best-fit effective 
radius R e , the Sersic index n, the ratio of minor to major axis b/a, and the position angle P.A. (in degrees East of North) 
correspond to the median of the distribution of results, and the uncertainties represent the 68% confidence intervals around 
the best fit. 



TABLE 5 

Best-fit Structural Parameters from the SINFONI Ha Maps and PSF-Matched NIC2 H [60 Band Images 



Object 


Map 


FWHM resolution" 


Re" 

(kpc) 




b/a b 


P.A. b 

(degrees) 


Q1623-BX528 


Ha 
H 160 


0"63 
0"63 


5.06+ 068 

4 Ql+l^f 

4 - yl -i 73 


29+ 69 
-Q 15 

u - ZJ -0 73 


53+ - 26 
" -0.10 
44 +H - 32 


_ 4 7 1+30.4 
HZ -'-6.5 
31 3+14.7 


Q1623-BX663 


Ha 
H 160 


Of'39 
0"39 


3.97+1-15 
3 ^7+8:^ 

J - J '-0.6S 


sss 

UJO -0 16 


72+°' 17 
0.82*1 


4.9+24.3 
"-34 


SSA22a-MD41 


Ha 

H 160 


o;'44 

0"44 


5 30+ " 
5 64+0-85 


0.28+o«o 

35+° :55 
-0.24 




33 2+ 9 ' 3 
36.7+];° 


Q2343-BX389 


Ha 
H 160 


Of'54 
0"54 


6.10 +L62 
5 99+o:l§ 

J -'^-0.90 


a 3 ,40.86 
u - z;, -0.13 


OAS"* 

47+0- 33 
u - -0.09 


-50.9+' 48 
-44 fi+'S:^ 


Q2343-BX610 


Ha 
#160 


Of'39 
0"39 


5 12+ 201 

4 07+0.62 


56+ 71 
u - ;,o -0.35 

45+ - 26 
u - HJ -0.20 


75+" 17 

66+ 22 
UOD -0.09 


11 9+20.1 
-14.9 

10 1+ 12 ' 5 


Q2346-BX482 


Ha 
#160 

Ha 

#160 


0"17 

o:'i5 

0!'50 

o;'5o 


6 07 +7 09 
D - u '-1.46 

6 22+°- 13 

601^3 

6 22 +2 u 
°- zz -0.91 


35 +3 - 47 

0.14±0.01 

27+ 70 
u - z '-0.19 

1 5+ 050 


59+ 21 

"""-0.12 

0.48 ±0.01 

64 +0 20 
50+ 017 


-50.2+ 32 , 4 

-60.9+|j 

~ 58 -U-ia2 

-55 4+ 13 - 7 
JJ - -8.0 



NOTE. — Structural parameters derived from single-component Sersic profile fits to the two-dimensional Ha maps from SINFONI 
and to the Hi^o band images from NIC2 convolved to the effective spatial resolution an d resa mpled to the pixel scale of the SINFONI 
data. The results are obtained from a series of 500 GALFIT runs as described in § 14.11 and the best fit value and uncertainties 
correspond to the median and 68% confidence intervals of the distribution of results in each parameter. 



a PSF FWHM of the SINFONI Ha maps, and of the convolved NIC2 H 160 band images. 

b Main structural parameters derived from the GALFIT fits: the effective radius radius i? c , the Sersic index n, the ratio of minor to 
major axis b/a, and the position angle P.A. (in degrees East of North). 



Kpc-scale rest-optical imaging of z ~ 2 galaxies 



TABLE 6 

Emission line contributions to the band flux densities 



Object / BB ([OIII]A5007) / BB ([OIH] A4959) / BB (H/3) / BB (tot) 



Q1623-BX528 

Q1623-BX663 < 0.05 < 0.02 <0.12 < 0.20 

SSA22a-MD41 0.11 ±0.01 0.043 ±0.005 0.023 ±0.004 0.17 ±0.01 

Q2343-BX389 0.19 ±0.02 0.050 ±0.008 < 0.054 < 0.29 

Q2343-BX610 0.034 ±0.004 0.013 ±0.002 0.018 ±0.002 0.064 ±0.005 

Q2346-BX482 0.062±0.006 0.019±0.003 < 0.037 <0.12 



NOTE. — Individual and total fractional contributions of the [O III] AA 4959, 5007 
and H/3 emission lines to the integrated broad-band flux density measured from the 
NIC2 #i6o band data. The uncertainties correspond to 1 a uncertainties from the line 
fluxes and broad-band magnitudes. For undetected lines, the 3 a upper limit on the 
line flux has been used. 
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BX528 
2=2.2683 




FWHM 



BX663 
2 = 2.4332 




FWHM 



MD41 
2=2.1704 




FWHM 



BX389 
2 = 2.1733 



FWHM 




BX610 
2 = 2.2103 




BX482 
2=2.2571 




FWHM 



Fig. 1. — HST NICMOS/NIC2 Hm maps of the galaxies, probing the broad-band emission around rest-frame 5000 A. The target name and 
Ha redshift are labeled in the top corner of each panel. The color coding scales linearly with flux density from white to black for the minimum 
to maximum levels displayed (varying for each galaxy). The FWHM of the effective PSF is indicated by the filled circle at the bottom left of 
each corner; the angular resolution is 0"145, or ~ 1.2 kpc at the median z = 2.2 of the sources. The angular scale of the images is shown with 
the l"-long vertical bar next to each galaxy. The black solid circle overlaid on each map shows the aperture used to extract the integrated Hieo 
band flux density (diameter of 3"). In all maps, North is up and East is to the left. 
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Fig. 2. — Broad-band optical to near-IR SEDs of the galaxies, including the Hieo flux density measured from the NIC2 data presented in this 
paper. The large white-filled circles and error bars represent the observed photometry and lcr uncertainties (Table \Q. The grey-filled circle 
shows the K band flux density corrected for the contribution by the Ha line emission, estimated from the available SINFONI data and given 
explicitly in each panel. The solid lines show the synthe tic spectrum of the bes t-fitting stellar p opu lation model, a ssuming a constant star 
formation rate, solar met allici tv. thelCalzetti et al.l ( 120001) reddening law, and the Chabrier (2003) or Kroupa (2001) IMF, obtained with the 
iBruzual & C harlot (2003) and Marastonl d2005l) models (blue and red lines, respectively). The best-fit age, Ay, stellar mass, and SFR for each 
galaxy and synthesis code are listed in each panel. 
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Fig. 3. — Stellar and dynamical properties of the galaxies observed with NIC2. The data points for the NIC2 targets are indicated with 
yellow, and the data for the comparison samples are taken from Forster Schreiber et al. (2009). (a) Stellar mass and star formation rate derived 
from SED modeling, compared to those of the SINS Ha sample at 1.3 < z < 2.6 and of A^-selected galaxies from the CDFS ( Wuvts et al. 
2008) in the same redshift interval and at A" s < 22.0 mag (the magnitude of the faintest of the SINS Ha sample galaxies in the K band). The 
SINS data points are shown with large filled dots, their projected distribution onto each axis with blue-hatched histograms, and their median 
magnitudes as blue-hatched horizontal bars. The CDFS data are plotted with small grey dots, grey filled histograms, and grey -hatched bars. The 
histograms are arbitrarily normalized. The galaxies classified as disk-like and merger-like by kinemetry ( Sh apiro et al. 120 08) are plotted as red- 
and green-filled circles. Sources that were known to host an AGN based on optical (rest-UV) or previous long-slit near-IR (rest-frame optical) 
spectroscopy are indicated with a 6-pointed skeletal star. SINS targets that were not detected in Ha line emission are marked as cyan-filled 
circles, (b) Half-light radius and circular velocity va derived from the SINFONI Ha observations, compared to those of the entire SINS Ha 
sample. The color scheme for data points is as for the left panel, and the projected distributions and median values for the SINS Ha sample are 
shown with the grey histograms and grey-hatched bars. The NIC2 sample spans roughly an order of magnitude in stellar mass and SFR, and 
lies at the high end of the distribution of Ha sizes and circular velocities inferred for the SINS Ha sample. 
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Fig. 4. — Results of parametric fits to the two-dimensional Hieo band surface brightness distribution of our NIC2 sample. The structural 
parameters (effective radius R e , Se rsic in dex n, axis ratio b/a, and position angle P.A.) were derived from a series of 500 single-component 
Sersic model fits, as described in § 14. 1.11 Each row shows the results for a given galaxy, with panels as follows. Far left: n versus /? c , where 
small grey dots show the results of the 500 iterations, histograms correspond to the projected distributions onto each axis, and the large black 
dot and error bars indicate the median and 68% confidence intervals adopted as best-fit parameters and 1 a uncertainties. Middle left: Same 
as the far left panel but for b/a versus P.A. Middle right: NIC2 images with overplotted ellipse of center, semi-major axis, axis ratio, and 
orientation equal to the adopted best-fit position, _/? e , b/a, and P.A. Far right: Residual map obtained by subtracting the NIC2 image with the 
Sersic model with best-fit parameters. All images are plotted with a linear intensity grey scale (with negative to positive residuals ranging from 
white to black). 
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Fig. 5. — Same as Figure[4]for the two-component Sersic fits for BX528 and BX389. In the three first panels to the left and for BX528, blue 
and red colors correspond to the results for the southeast and northwest components, respectively. For BX 389, blue is used for the main part of 
the galaxy, and red for the small southern companion. 
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Fig. 6. — Radial profiles in Him band light of our NIC2 sample. The profiles are azimutally-averaged along concentric ellipses with center, 
semi- majo r axis, axis ratio, and orientation equal to the best-fit position, R e , b/a, and P.A. values from the Sersic model fits to the galaxies 
(see § |4,lt . Each row shows the results of a given galaxy, with the panels as follows. (Left) Mean intensity in units of [ijy as a function of 
semi-major axis radius of the elliptical annuli; the black curve and dots show the profile extracted from galaxy's image while the orange curve 
shows the corresponding profile for the best-fit Sersic model. The radial intensity profile of the PSF (along concentric circular annuli) is plotted 
in green for reference. (Middle) Same as for the left panel for the mean surface brightness in units of mag arcsec~ 2 , with the best-fit structural 
parameters labeled explicitly. (Right) Ratio of the mean intensity profiles of the data and of the best-fit model as a function of semi-major axis 
radius of the elliptical annuli. For all galaxies except BX 528, single-component Sersic models are considered, and the derived effective radius 
is indicated by a vertical dashed line. For BX 528, the model curves plotted here are extracted from the image for the best-fit two-component 
model but along elliptical annuli defined by the best-fit single-component model. The different intensity units between the left and middle 
panels are chosen to highlight in turn the central brighter and outer fainter regions of the sources. 
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Fig. 7. — Morphological parameters for SINS galaxies derived from the NIC2 Hu,o band images. Each galaxy was analyzed as if it was placed 
at the maximum redshift of the sample, i.e., z = 2.43. (Left) Gini versus ^ (Multiplicity). (Right) Gini versus M20. The AO-assisted SINFONI 
Ha map of BX482 was also analyzed in the same manner, and the results are plotted (star symbol) along with the NIC2 results (squares) in 
both panels. 



MD41 

2 = 2.1704 





Fig. 8. — Comparison of rest-frame UV and optical morphologies for MD41. North is up and east to the left in both images, and a bar of 1" 
in length indicates the scale in the images. (Left) ACS !si4 band map at the original angular resolution of the data, with the FWHM of the PSF 
indicated by the filled circle at the bottom left corner. (Middle) Same as the left panel for the NIC2 Hiea band map. (Right) Color-composite 
after accurate PSF matching and registration of the ACS and NIC2 maps; in this RGB image, the red channel is assigned to the Hm® band, and 
the green and blue channels to the f'gi4 band. The colors for all images plotted correspond to a linear flux density scale. The white cross shows 
the location of the geometrical center of MD 41 as determined from the NIC2 map. Overall, the images are strikingly similar, although on small 
scales there are noticeable differences. In particular, the two brightest clumps in (gi4 band on the southwest edge of the galaxy are much fainter 
in Hmx, band and, consequently, have bluer colors than the bulk of the source. In addition, two compact regions appear in the (814 band image 
alone, which are offset from the main body of the galaxy. The colors of these compact offset sources are thus significantly bluer than that of the 
main part of the galaxy. These sources also have no counterpart in line emission within ±1000 km s -1 of the Ha line of MD 41 in our SINFONI 
data. Therefore, we cannot determine whether they are physically associated with MD41. 
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Fig. 9. — Comparison of fligo band images and Ha line maps and velocity fields from SINFONI for three of our NIC2 targets: BX528, 
BX663, and MD41. Each row corresponds to a given galaxy, with the panels as follows. (Left) flieo band map at the original NIC2 angular 
resolution. (Middle) Ha line map with contours of the PSF-matched NIC2 H\$a band image overplotted in light grey. (Right) Ha velocity field 
with contours of the PSF-matched NIC2 H\(,o band image overplotted in white; the velocities are color-coded from blue to red for blueshifted 
to redshifted line emission relative to the systemic velocity and over the range labeled at the top left. In each panel, the FWHM of the PSF of 
the corresponding image is shown by the filled circle in the bottom left corner. The vertical bar indicates an angular scale of 1" (~ 8 kpc at the 
redshift of the sources). The white crosses mark the geometrical center of the sources. The total on-source integration times are also given for 
the NIC2 and SINFONI data, as well as whether the SINFONI data were obtained without AO or with AO using a natural or the laser guide star 
("NGS-AO" and "LGS-AO," respectively). 
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Fig. 10.— Same as Figure^for the three other NIC2 targets: BX389,BX610, andBX482. 
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Fig. 11. — Comparison of the structural parameters from single-component Sersic profile fits to the two-dimensional H\(q band and Ha 
surface brightness distributions. For consistency, the fits were performed on the Hi&> band images matched to the PSF and resampled to the 
0"125 pixeP 1 scale of the lower ~ 0"5 resolution Ha maps. Data for the different targets are color-coded as indicated in the leftmost panel, (a) 
Effective radius, (b) Sersic index, (c) Axis ratio, (d) Position angle (in degrees East of North). The solid line indicates the one-to-one relation. 
The parameters derived from the Ha and smoothed Him band maps at w 0"5 resolution are in good agreement and typically within their 1 a 
uncertainties. They also agree well with the values obtained from fits to the original higher resolution NIC2 images (see Table|5j- 
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Fig. 12. — Radial light profiles of the Ha and H[6o band emission of the targets. The profiles and best-fit single-component Sersic model 
overplotted are based on the ~ 0"5 resolution SINFONI Ha maps (using the 0"l25 pixel -1 scale, taken under seeing-limited conditions or 
with NGS-AO for BX663) and PSF-matched NIC2 H i60 images except for BX482, where the AO-assisted SINFONI and original resolution 
NIC2 data both at 0."05 pixel -1 are used (see §|5>. Fr om left to right, each row shows for a given source the mean intensity in lin ear flux units, 
the surface brightness in magnitude units, and the ratio of observed intensity profile to that of the best-fit Sersic model (§ 15. 2t . The data are 
plotted with solid lines and dots, and the models are shown with dashed lines; blue and red colors are used for Ha and Hk,o band, respectively. 
The profiles are computed along ellipses with center, axis ratio, and P. A. corresponding to the Sersic model that best fits the NIC2 images; 
the maximum radius is determined by the effective field of view of the SINFONI data and so as to avoid the noisiest parts at the egdes. The 
intensity and surface brightness profiles are all normalized to unity and 25 mag, respectively, at the central position. The ratio of observed and 
model intensity profiles are as computed before normalization, and the horizontal grey line in the right-hand panels indicates a ratio of unity. 
The vertical dotted lines in all panels show the best-fit effective radii from the Sersic fits, and the SINFONI PSF radial intensity profiles are 
plotted in the left-hand panels in green for reference. The close correspondence between Ha and rest-frame ~ 5000 A radial light distributions 
is striking, reflecting the similarity in the maps shown in Figuresl9landllOI 
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Fig. 13. — Sersic index n versus effective radius i? c for our six SINS galaxies and comparison samples at similar redshifts but with different 
selection criteria. All galaxies were observed with the same NICMOS NIC2 camera through the same F160W filter and the Sersic model fits 
were performed in the same w ay for all samples, allow ing consistent comparisons, (a) n versus R c fo r the SINS and the A"-selected samples 
(labeled "vD" for objects from Ivan Dok kum et al. 2008 and "Kriek" for those from Kriek et al. 2009). The images used for all sources were 
drizzled to a 0."038 pixel scale. Results for the SINS galaxies based on the full 4-orbit data and 3-orbit subsets (as 17 of 19 of the AT-selected 
objects) are indistinguishable, so the parameters derived from the 4-orbit data are plotted. The SINS targets are plotted with filled circles 
following the color scheme as gi ven by the labels . Squares and triangles correspond to the quiescent and star-forming ^-selected samples of 
Ivan Dokkum et alH d2008t) and Kriek et al. (2009), respecti vely. The A GN sources are distinguished with a starred circle for BX663 in our 
SINS sample, and with grey- filled triangles for those of the Kriek etaD sample, (b) Same as previous panel for the SINS and 24 /im-selected 
dusty IR-luminous sample of Da syra et al. I (2008; labeled "Dasyra"). All images used here were drizz led to 0"076 pixels. Large and small 
symbols are used for results based on 2- and 1-orbit integrations, respectively. Ten of 11 sources from lDasvra et al.l f2008) show signatures 
of AGN activity in their mid-IR spectrum, either dominant or along with PAH features indicative of star formation; the only star-formation 
dominated object is the second largest with second highest n (4.06 kpc and 2.59). 
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Fig. 14. — Morphological parameters derived from NIC2 Hi6o band imaging for our SINS targets and comparison samples at similar redshifts 
but with different selection criteria. Each galaxy was analyzed as if it was placed at the maximum r edshift of the combined samples, i.e. 
z = 2.60. (Top, Left) Gini versus ^ for the SINS and /("-selected samples (labeled "vD" for objects from van Dok kum et al.l l2008 and "Kriek" 
for objects from Kriek et al. 2009). As in Figure[T3] AT-selected objects are indicated with solid squares (vD) and triangles (Kriek). For vD 
objects, medium-sized symbols are used for 3-orbit exposures and small symbols for 2-orbit exposures, while, for Kriek objects, all points 
represent 3-orbit exposures. SINS objects are indicated with open squares, with the largest symbols for 4-orbit exposures, medium-sized 
symbo ls for 3-orbit subse ts of the data, an d small symbols for 2-orbit subsets. Our NIC2 images were drizzled to the 0''038 pixel sampling 
of the Ivan Dokkum et akl and Kriek et al. data for consistent comparison. The non-parametric statistics for the empirical PSF constructed 
from stars in our 0."038 pixeF 1 NIC2 images are indicated with an open circle. (Top, Righ t) Gini versus M20 for the SINS, vD, and Kriek 
samples. Sym bols are as in the previous panel. (Bottom, Left) Gini versus ^ for the SINS and D asvra et al.l (2008; labeled "Dasyra") samples. 
Dasvr a et al. I objects are indicated with solid black squares, with medium-sized symbols for 2-orbit exposures and smaller symbols for 1-orbit 
exposures. Again, SINS objects are indicated with open squares, with the largest symbols for 4-orbit exposures, medium-sized symbols for 
2-orbit subse ts of the data, and smallest symbols for 1-orbit subsets. Our NIC2 images were here drizzled to the 0"076 pixel scale of the 
iDasvra et al.l data, and the non-para metric statistics for t he empirical PSF from these images are also indicated with an open circle. (Bottom, 
Right) Gini versus M20 for SINS and Dasv ra et al.l 0008) samples. Symbols are as in the previous panel. 
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Fig. 15. — Struct ural parameters of our six SINS galaxies and of the ^-selected comparison sample from van Dokkum et al. (2008, "vD") and 
iKriek et al.l < f2009l "Kriek"). plotted as a function of global stellar properties, (a) Effective radius i? c versus stellar mass M». (b) Sersic index n 
versus M„. (c) W c versus specific star formation rate s SFR. (dj n versus sSFR. Symbols and colors are the same as in Figure[T3] and labeled in 
panel (a). All stellar masses are for a Chabrier (2003) IMF. In the range of M* where our SINS objects and the A"-selected samples overlap, the 
data reflect the wide range in size (and in n) that is known to exist for galaxies of similar stellar masses and out to z ~ 3. Size and star formation 
activity are clearly correlated, again as seen for the massive galaxy population locally and at high redshift. 
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Fig. 16. — Same as Figure[T5]for the non-parametric morphological coefficients. The top row shows morphological coefficients plotted as a 
function of M„, while the bottom row shows them as a function of specific SFR. As for the structural parameters (R e and n), a large diversity 
of non-parametric morphological coefficients are found in the region of overlapping M t , while significant correlations are found between 
morphology and specific SFR. 
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Fig. 17. — Same as Figure [15] for the specific star formation rate as a function of the stellar mass surface density. Half of the total stellar 
mass is assumed to be enclosed within a radius corresponding to the circularized effective radius /fecta = yb/a, givin g estimates of the 
projec ted s urface densi t y. Th ere is a clear trend outlined by our SINS targets and the A"-selected comparison samples of van Dokku m et al.l 
(2008) and Krie k et al.l ( I2009I) . which span overall three orders of magnitude in stellar mass surface density and even more in specific SFR (a 
similar trend is seen when using the R c instead for estimates of the deprojected sur face density). For comparison with the more general massive 
galaxy population at z ~ 2, the median relationship derived bv lFranx et al.l (f2008) from a large AT-selected sample in GOODS-S and based on 
ground-based near-IR imaging to derive rest-frame optical sizes is overplotted as a grey thick line. Our SINS and the A"-selected comparison 
samples consistently follow this relationship. 
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Fig. 18. — Color, observed mass-to-light ratio, and stellar mass maps of MD41. North is up and east to the left in all panels, and a bar of 
1" in length (approximately 8 kpc at the redshift of the source) indicates the scale in the images. (Left) Observed ;si4-//i60 colors from the 
ACS/F814W and NIC2/F160W images, matched to the same PSF with FWHM indicated by the filled circle at the bottom left corner. (Middle) 
Observed stellar mass to rest-frame g-band light ratio, computed from the 1^14- Hu,o colors using the relationship derived in Appendix A. (Right) 
Stellar mass surface density, derived from the M*/Lg Kt ratio of the previous panel and the NIC2 Huo band map. Pixels with a S/N > 3 in each 
of the ACS/F814W and NIC2/F160W images are shown, while others are masked out. 
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Fig. 19. — Relationship between observed igu-Him colors and the ratio of stellar mass to rest-frame g-band luminosity M»/Lj, csC . Both the 
colors and the M ±/Lf st correspond t o obse rved quantities, uncorrected for t he effects of dus t extinction. T he various curves plo tted in colors are 
computed from Bruzual & Chariot ( 2003) models with solar metallicity, a Chabrier (2003) IMF, and the Calzetti et al. ( 2000) reddening law. 
Different colors are used for different star formation histories: red to blue represent SFHs of increasing star formation timescale from a single 
stellar population, a suite of exponentially declining SFRs with e-folding timescales between 10 Myr and 1 Gyr, and constant SFR. Different 
linestyles are also used for models computed with different values for extinction: Ay = (solid), 1 (dashed), 2 (dotted), and 3 mag (dash-dotted). 
Age increases along each model curve from blue to red hu— H\$a colors and low to high M*/Lj cst ratios. As shown by the arrow, the effects 
of extinction are roughly parallel to the locus of the various models. The colors are synthesized from the model spectrum redshifted to the 
z = 2.1704 of MD41, using the transmission curves for the ACS F814W and NIC2 F160W filters. The large white-filled circles and thick black 
line shows the mean relationship derived from the solar metallicity models. For comparison, the smaller grey-filled circles and thinner black 
line is the relationship that is derived from models with 1/5 solar metallicity. The grey-shaded histogram shows the distribution of observed 
i8i4 -#160 colors of all individual pixels with S /N > 3 in each band for MD 41 (see § 17. It . 



